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^ Context. Effective temperatures of early-type supergiants are important to test stellar atmosphere- and internal structure-models of massive and 
, intermediate mass objects at different evolutionary phases. However, these T^if values are more or less discrepant depending on the method 
' used to determine them. 

, Aims. We aim to obtain a new calibration of the T^ff parameter for early-type supergiants as a function of observational quantities that are: a) 
^ highly sensitive to the ionization balance in the photosphere and its gas pressure; b) independent of the interstellar extinction; c) as much as 
possible model-independent. 

Methods. The observational quantities that best address our aims are the (Ai,D) parameters of the BCD spectrophotometric system. They 
T-H ■ describe the energy distribution around the Balmer discontinuity, which is highly sensitive to T^ff and logg. We perform a calibration of the 
lO ' (A-uD) parameters into T^ff using effective temperatures derived with the bolometric-flux method for 217 program stars, whose individual 

uncertainties are on average lAreffl/T^u = 0.05. 

Results. We obtain a new and homogeneous calibration of the BCD (Ai,D) parameters for OB supergiants and revisit the current calibration 
of the (A[,D) zone occupied by dwarfs and giants. The final comparison of calculated with obtained Teff values in the (Ai,D) calibration show 
that the latter have total uncertainties, which on average are ^T^^/T'^g- ±0.05 for all spectral types and luminosity classes. 
ILJ Conclusions. The effective temperatures of OB supergiants derived in this work agree on average within some 2 000 K with other determinations 
. ^ found in the literature, except those issued from wind-free non-LTE plane-parallel models of stellar atmospheres, which produce effective 
temperatures that can be overestimated by up to more than 5 000 K near T^g = 25 000 K. 

Since the stellar spectra needed to obtain the (Ai,D) parameters are of low resolution, a calibration based on the BCD system is useful to study 
stars and stellar systems like open clusters, associations or stars in galaxies observed with multi-object spectrographs and/or spectro-imaging 
devices. 
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1. Introduction 

The effective temperatures of early-type stars of luminos- 
ity classes V to 111 are similar whatever the method used 
to determine them. Among the most commonly used meth- 
ods are those based on: fits of the observed absolute spectral 
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energy distributions (hereafter ASEDs) with iKurucz d 19791) 
line-blanketed stellar atmosphere models; line profile fit- 
tings; Stromgren and Geneva photometr ic color indice s. Using 
Breger's spectrophotometri c catalogue ('BregeiJ' l976al) for the 
ASED in the visual range, Morossi & Malagnitijl 1985 ) and 
Malagnini & Morossil d 19901) obtained Tgs values with inter- 



nal uncertainties below 5%. The uncertainties are of the or- 
der of 10%, either when the TDl far-UV fluxes and lUE low 
resolution spectra a re used in combination with fluxes from 
Breger's catalogue (iMalagnini et al J 1 1 983 . 119861; iGulati et alJ 
1989), or when H and He absorption l ine pr ofiles are fit- 



ted with models dMorossi & Crivellari Il980l) . In general. 



Teff values derived using calibrated Stromgren and Geneva 
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photometric indices present low uncertainties (Balonal 1984 : 
Moon & Dworetskvlll985l:ICasteUilll99ll:lAchmad et al.lll993b . 



However iNapiwotzki et alj (11993 ^, based on a critical com- 
parison of several calibrations of Stromgren intermediate -band 
uvby-y8 photo metric indices, recommende d the use of the cali- 
bration done by lMoon & Dworetskyl(ll985h . corrected for grav- 
ity deviations. 



As regard s B supergiants, ICode et al. (Il976h and 



Underbill et al. (1 19791) derived effective temperatures based 



on spectrophotometric observations in the far-UV, visible and 
near-IR spectral regions. Later, effective temperature determi- 
nations for these stars were made using the silicon lines in 
the optical spectral region, either by fitti ng the line profiles or 
by measuring their equivalent widths (iBecker & Butlei|[l990l 
McErlean et al] Il999t iTrundle et all |2004 . Nowadays, meth- 
ods based on adjustments of He i and He n fine p rofiles, and/or 
the line intensity ratio Si l y/Si iii are preferred (IHerrero et al. 



,2002'; 'Repolust et al."2004l, 'Martins et al."2005"; 'Crowther et al 
2006; Benaglia et al. 2007; Markova & Puis 2008; Seade et al. 
20081) . Recent estimates of the effective temperature of early 
B-type supergiants have shown that the values obtained with 
non-LTE blanketed models including winds (hereafter non- 
LTE BW models) are systematically lower than those de- 
rived with models w i thout winds (he reafter wind-free models) 
jLefever et al.l 120071: ICrowther et al. I [2006). Differences range 
roughly from to 6 000 K and they tend to be lower the later 
the B -sub-spectral type ("Marko va & PulslbOOSl) . On the other 
hand, Morossi & Crivel lari (198C^ noted that the effective tem- 
peratures of B-supergiants derived with photometric methods 
are in general higher than those obtained spectroscopically. 
Physical characteristics and phenomena like activity and/or in- 
stabilities taking place in the extended atmospheric layers af- 
fect more significantly the spectral lines than the continuum 
spectrum. Thus, the genuine signature due to the Teff caiTied 
by the lines could be somewhat bluiTed. 

The calibrations of effective temperatures for B supergiants 
are of great importance since these stars are in a significant 
phase of the evolutionary sequence of massive stars. They are 
also the main contributors to the chemical and dynamical evo- 
lution of galaxies. Accurate effective temperatures are then 
needed to construct HR diagrams and to test the theories of 
stellar structure and evolution, as well as to estimate the chem- 
ical content of the stellar environment. Effective temperatures 
are also necessary to study the physical processes in the atmo- 
sphere, such as non-radial pulsations or stellar winds (radia- 
tive forces; changes in ionization). Regarding the stellar winds, 
the effective temperatures are particularly useful in discussing 
terminal velocities, mass loss rates, the bi-stability jum p, and 



the \yind momenturn luminosity relationship (.K udritzk i et al 
2003HCrowther et al.ll2006HMarkova & Pulsll2008l) . 



As a consequence of the large discrepancies found in 
the Teff estimates of B supergiants, the current temperature 



scale is being revisited jMarkova & Pulsl l2008t ISearle et al 



20081) . In this context, we present an independent and homo 
geneous temperature calibration for B-type dwarfs to super 
gian ts, based on the use of the BCD spectrophotometri c sys 



based on : a) measurable quantities that are strongly sensitive 
to the ionization balance in the stellar atmosphere and to its gas 
pressure, thus being excellent indicators of Teff and log g; b) pa- 
rameters that describe the visible continuum spectrum, whose 
atmospheric formation layers are on average deeper than those 
for spectral lines. 

Our first step will be to determine the effective temperatures 
of our sample of Galactic B-type supergiants. However, to per- 
form a consistent calibration of the BCD {Ai,D) parameters for 
supergiants, we re-determine the calibration into effective tem- 
perature of the BCD {A\,D) domain coiTesponding to B-type 
dwarfs and giants (D = size of the Balmer jump; Ai = mean 
spectral position of the Balmer discontinuity (BD); see further 
explanations on these parameters in ^and Appendix jAt . We 
use a large and homogeneous sample of B stars observed in the 
BCD system, and newly-derived effective temperatures for all 
of them, based on the bolometric-flux method (hereafter BFM). 
We leave for another contribution the discussion of BCD caU- 
brations related to log g, visual and bolometric absolute magni- 
tudes. 

The present paper is organized as follows: In ^we briefly 
describe the BCD spectrophotometric system and the advan- 
tages of its use. In ^we present the BFM on which the deter- 
minations of the stellar effective temperature and angular diam- 
eter (6) are based. Observations and the Tetf values determined 
with the BFM are presented in ^ The uncertainties of the ef- 
fective temperatures and angular diameters obtained with the 
BFM are discussed in ^ Comparisons of our T^ff and 6 deter- 
minations with those obtained by other authors are given in |6] 
In ^ we present the empirical temperature calibration curves 
and discuss the accuracy of the Tf.ff(Ai,D) values obtained. A 
discussion and global conclusions are presented in ^ and ^ 
respectively. 

2. The BCD system 

The Paris spectrophotometric classification system of stellar 
spectra, b est known as BCD (Barbie r-Ch alonge-Divan), was 
defined by iBarbier & Chalongd (1194 ih and IChalonge & Divan 
(1952t). The original presentation of this syst em is given i n 
French; explanat i ons in English can be f ound i n Dufav (1964 ); 

ibivan, (,19921) . 



Underhilll (Il966l) ; lUnderhill & DoazanI (Il982b ; 



tem ( iBarbier & Chalongd 194"lt IChalonge & Divanll952l) . This 
method has numerous advantages (see mainly because it is 



A short overview of the system is given in Appendix ^A] The 
BCD system is based on four measurable quantities in the con- 
tinuum spectrum around the BD: D, the Balmer jump given in 
dex and determined at A31QQ A; Ai, the mean spectral posi- 
tion of the BD, usually given as the difference A[-31QQ A; Oj,,., 
the gradient of the Balmer energy distribution in the near-UV 
from /13100 to A31QQ A, given in yum; the gradient of the 
Paschen energy distribution in the wavelength interval AA4QQQ- 
6200 A given in fim. The sole BCD parameters that are relevant 
to the present work are: D, which is a strong function of Teff, 
and Al that is very sensitive to log g. 

The use of the {Ai,D) pair to determine the spectral clas- 
sification and the stellar fundamental parameters presents nu- 
merous advantages, not only because the BD is a well visible 
spectral characteristic for stars ranging from early O to late F 
spectral types but also because: 
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a) The parameters (Ai ,D) are obtained from direct measure- 
ment on the stellar continuum energy distribution. This implies 
that, on average, they are relevant to the physical properties of 
photospheric layers which are deeper than those described by 
spectral lines; 

b) Each MK (Morgan & Keenan) spectral type-luminosity 
class (SpT/LC) is represented by wide intervals of A\ and D 
values, which implies high SpT/LC classification resolution: 
Al ranges from about 75 A for dwarfs to -5 A for supergiants, 
while D ranges from near 0.0 dex, for the hottest O stars and 
F9 stars, to about 0.5 dex, for the A3-4 stars; 

c) Typical Icr measurement uncertainties affecting D and 
Al are: 6D < 0.02 dex and 6A\ < 2 A, respectively. Thus, from 
b) and c) we find that hardly any other classification system 
has reached such a high resolution, especially concerning the 
luminosity class; 

d) The parameter Ai is independent of the interstellar 
medium (ISM) extinction, while D has a low E(B-V) color 
excess dependence, roughly 6D - 0.03 E(B-V) dex, which 
is almost insensitive to the selective absorption ratio Ry - 
Av/E(B - V). The 6D difference is produced by extrapolation 
of the Paschen energy distribution from /14000 A to A3700 A, 
which caiTies the ISM reddening of the Paschen continuum. 
The low ISM extinction dependence is however of great inter- 
est for the study of early-type stars, since they are frequently 
distant and strongly reddened; 

e) Spectra needed to obtain the iAi,D) measurements are of 
low resolution (8 A at the BD). This means that numerous faint 
stars can be observed with short exposure times. Furthermore, 
the flux calibrations required to derive the BCD parameters rely 
on spectral reduction techniques which are easy-to-use and of 
common practice; 

f) The BCD system is generally used for 'normal' stars, 
i.e. objects whose atmospheres can be modeled in the frame- 
work of hydrostatic and radiative equilibrium approximations. 
However, since both the photospheric and the circumstellar 
components of the BD are spectroscopically well separated, it 



stars dDivan & Zored 1982; Zorec 


198(4 Zorec &Briotl 1991: 


ChauviUe et al. I2OOII; Zorec et alJ 


2OO5I 


Vinicius et al.j2006b; 



chemically peculiar stars (He-W group) tCidale et al.ii2OO70 . 

3. The bolometric-flux method (BFM) 

To obtain the calibration of the BCD iAi,D) parameters into ef- 
fective temperature, we could simply adopt for each star the av- 
erage of all the Tgff values found in the literature. Nevertheless, 
these quantities were obtained with several heterogeneous 
methods which carry more or less systematic differences on 
the estimates of Teff. Since most T^ff determinations are based 
on adjustment of the observed ASEDs with theoretical ASEDs, 
or on fittings of line profiles with models, the properties of the 
studied stars are implicitly assumed to be in accordance with 
the physical characteristics of the best fitted stellar model at- 
mosphere. Instead, in the present contribution, we preferred to 
determine the effective temperature by using the total amount 
of radiated energy, so that the details of its distribution are of 



marginal importance and the dependence on models of stellar 
atmospheres are kept to a minimum. While models of stellar 
atmospheres produce Balmer jumps which are close to the ob- 
served ones, the theoretical Ai parameter may differ somewhat. 
In fact, Al depends on the distribution of the emergent radiation 
fluxes near the limit of the Balmer line series, where the theo- 
retical uncertainties concern the treatment of the non-ideal ef- 



fects in the hydrogen upper level populations (iRohrmann et al 
2003h. 



We decided to adopt a single method for all program stars 
and to estimate their effective temperatures using its definition, 
where the incidence of the model-dependence is in principle 
strongly minimized. 

By definition, the effective temperature of a star is: 



' eff 



(1) 



where ctr is the Stefan-Boltzmann constant; / is the stellar 
bolometric radiation flux received at the Earth, corrected for 
the ISM extinction; is the angular diameter of the star If 
the radiation field coming from the stellar interior were the 
sole energy source in the atmosphere, the effective tempera- 
ture deduced from the bolometric flux would be the same as 
that derived from the analysis of stellar spectra with model 
atmospheres in radiative and hydrostatic equilibrium. If the 
bolometric flux / and the angular diameter 6 were issued en- 
tirely from observations, the T^ff deduced from Eq. ([TJ could 
then also be considered a genuin e observational par ameter. The 
BFM was applied in this way by Code et al. ( 19761) using stel- 
lar fluxes observed with the OAO-2 satel lite and the angular di- 



amete rs determined interferometrically (IHanburv Brown et al 



19741) . Unfortunately, this could not happen in our case, as we 
do not have observed fluxes over the entire spectrum, nor do we 
have angular diameters measured for all program stars. In what 
follows the effective temperature and the angular diameter de- 
rived with the 'bolometric flux method' are called T^g and 6^, 
res pectively. 



Blackwell & ShallisI ( 119771) showed that a stellar angular 



diameter can be well reproduced with observed fluxes in near- 
IR and model atmospheres. The monochromatic stellar angular 
diameter is thus given by: 



(2) 



where /° is the absolute monochromatic flux received at the 
Earth corrected for the ISM extinction, - nFji is the emit- 
ted monochromatic flux at the stellar surface, is the so called 
'astrophysical' flux predicted by a model atmosphere. To repre- 
sent F,i we ha ve used the grids of ATLA S9 model atmospheres 
calculated by Castelli & Kuruczl ( 2003 ). In wavelengths lying 
in the Rayleigh-Jeans tail of the energy distribution, not only 
are the theoretical fluxes nearly independent of model charac- 
teristics, but the observed ASEDs also are mildly affected by 
the ISM extinction. However, they have the inconvenience of 
being frequently maiTed by infrared flux excesses of non stel- 
lar origin. We decided to calculate Eq. ^ in the red extreme of 
the Paschen continuum: AA 0.58-0.8yum. In this spectral region. 
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the theoretical radiative fluxes are still only sUghtly dependent 
on the particular model characteristics and the layers where this 
radiation field is formed have local electron temperatures close 
to Tgff, so that the color temperature of the energy distribution 
of OB stars in this wavelength interval also approaches Teff 
closely. While Tefj derived with Eq. ([TJ is insensitive to log g 
within the characteristic uncertainties of its determination, the 
parameter 6,i does depend slightly on models. Therefore, in 
Eq.© model atmospheres giving are chosen for gravity 
parameters log g - log g(J3, Teff), where /3 is the Hy6-line index 
of the uvby -;6 Stromgren photometry. Th e /3 parameters used 
are from the iHauck & Mermilliod ( 1975') compilation and the 
log giJ3, Teff) relations used are from[Castelli & Kurucz (2006). 
An ex tensive use of the BFM was made by lUnderhill et al] 
19791) . where unfortunately Eqs. ([TJ and ^ were iterated only 



twice, which left the derived eff'ective temperatures strongly 
correlated with their initial approximate values. 

3.1. Effective temperatures of B-type dwarfs to giants 

The ASEDs of our program stars, taken from the literature to 
calculate /, are in most cases observed in the wavelength in- 
terval ranging from Aa ~ 1200-1300 A to some At, in the far- 
IR. Many times they present intermediate gaps, according to 
the observational method or the instruments used. On the other 
hand, the angular diameter has been measured only for a few 
bright stars. Nevertheless, as in Eq. ([T]) Tgff depends on /'^'*, 
the effective temperature can still be reliably determined even 
though the unobserved spectral regions are represented using 
'modestly realistic' model atmospheres, as discussed in ^ The 
calculation of the effective temperature for dwarfs to giants is 
based on the following iteration: 

1) Adopt initial values of T^^. and log g; 

2) Interpolate the model ffuxes in the far-UV and IR spec- 
tral regions for the adopted (T^ff, logg) values ; 

3) Calculate the angular diameter 6^ using relation (|2]l in the 
near IR spectral region as detailed in the explanation following 
relation (|2]l, and assume 9-^ independent of the wavelength; 

4) Calculate the bolometric ffux / as follows: 



/ 

/obs 

6 



= /obs x[l+5] 
- dvv + Sir 



(3) 



JlR' 



UV 



where /obs and 6 represent, respectively, the contribution of the 
observed and unobserved spectral regions (extreme-UV, far- 
UV and IR) to the bolometric flux; Aa and Ah delimit the spec- 
tral range of ASEDs actually observed; /" and have the 
same meaning as given in Eqs. ([T]i and 

5) Introduce the bolometric flux given by (|3j into relation 
([TJ to obtain a new estimate of T^^ and accordingly, of log g - 
\ogg{Ti^,p)- 

6) Use the new estimates of {T^^^,\ogg) to continue the it- 
eration in step 2). 



0.4 



0.3 



0.2 



0.1 

o 
o 

+ 0.0 

COS 

O -0.1 



-0.2 



-0.3 



-0.4 



r 



HD41117 
• TD-1 data 

E(B-V)=0.494 mag 



Model: Tef/ = 20740,log g=2.5,7=1.00 
Model: Teff(7)=21100,log g=2.5,7=1.53 



1400 1600 1800 2000 2200 2400 2600 
A(A) 

Fig. 1. Fitting of the far-UV observed energy distribution of HD 
41117 with y-unmodified model fluxes for T^^^ (pointed line), 
and with y-modified model for Tf.ff(y) (dashed line). 



The iterations were performed until the difference between 
two consecutive T^^ values was smaller than 1 K. Depending 
on the star, this implies roughly 10 to 30 iterations. From ([3j it 
is obvious that the estimates of the effective temperatures may 
in principle be more uncertain the higher the value of 6; i.e. 
for the hottest stars. In Table[T]are listed the fractions <5uv and 
Sir for different effective temperatures and gravities. As seen 
in this table, S has a low dependence on log g. 



3.2. Effective temperatures of B-type supergiants 

The theoretical ASEDs predicted for supergiants by the wind- 
free plane-parallel model atmospheres with the effective tem- 
peratures issued directly from Eq. ([TJ, do not always fit well to 
the observed ASEDs in the near- and far-UV. Generally, the 
predicted fluxes are higher than the observed ones in wave- 
lengths A :$ 2200 A, while they are lower in the near-UV. For the 
obtained Teff they yield somewhat larger Balmer discontinuities 
than those observed. This might be partially due to the plane- 
parallel approximation of model atmospheres that probably is 
not suited to the extended atmospheric layers of these stars, to 
an incomplete treatment of the spectral line formation in such 
diluted atmospheres, to the omission of the e ffects produced on 
the photosphere by the stellar winds ( Abbott & Hummer 19^51: 



Gabler et al.ll989HSmith et al.l2002l:lMorisset et al.l2004 . and 

also to an insufl icient line blocking in th e model atmospheres, 
as discussed bv iRemie & LamersI ( Il98 1). Thus, to obtain the 
parameter of supergiants, we have slightly modified the use 
of relations ([TJ and ((2j through the following iteration proce- 
dure: 

1) Adopt approximate values of Teff and logg; 
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Table 1. Far-UV and IR fractions of the unobserved bolometric flux as a function of Teff and logg, the function A(Tgf[,logg), 
and temperatures as a function of Tgff(y). 



logg = 4.0 logg = 3.0 logg = 2.5 



reff(r) 






A[7'eff (r),g] 


r 


eff 






A[7'eft(r),g] 


r 


cff 






A[7'eff (r),g] 


r 


cff 


10000 


0.023 


0.059 


0.024 


1.5 
2.0 


9970 
9940 


0.024 


0.059 


0.028 


1.5 
2.0 


9970 
9930 


0.023 


0.059 


0.029 


1.5 
2.0 


9960 
9930 


12500 


0.105 


0.034 


0.088 


1.5 
2.0 


12360 
12220 


0.099 


0.034 


0.086 


1.5 

2.0 


12360 
12220 


0.093 


0.034 


0.081 


1.5 
2.0 


12370 
12240 


15000 


0.230 


0.022 


0.152 


1.5 
2.0 


14710 
14390 


0.218 


0.022 


0.135 


1.5 
2.0 


14740 
14470 


0.202 


0.022 


0.122 


1.5 

2.0 


14770 
14520 


17500 


0.386 


0.016 


0.182 


1.5 
2.0 


17090 
16640 


0.353 


0.016 


0.155 


1.5 
2.0 


17150 
16780 


0.318 


0.016 


0.134 


1.5 
2.0 


17200 
16880 


20000 


0.552 


0.012 


0.194 


1.5 
2.0 


19450 
18950 


0.484 


0.012 


0.157 


1.5 
2.0 


19600 
19160 


0.450 


0.012 


0.135 


1.5 
2.0 


19650 
19290 


22500 


0.706 


0.009 


0.191 


1.5 
2.0 


21940 
21340 


0.596 


0.010 


0.139 


1.5 
2.0 


22090 
21670 


0.546 


0.010 


0.109 


1.5 
2.0 


22190 
21860 


25000 


0.843 


0.008 


0.178 


1.5 
2.0 


24420 
23800 


0.693 


0.008 


0.106 


1.5 
2.0 


24670 
24310 


0.635 


0.008 


0.071 


1.5 
2.0 


24780 
24540 


27500 


0.955 


0.007 


0.153 


1.5 
2.0 


26950 
26380 


0.823 


0.007 


0.093 


1.5 
2.0 


27180 
26840 












30000 


1.060 


0.006 


0.126 


1.5 
2.0 


29520 
29010 


1.023 


0.006 


0.087 


1.5 
2.0 


29670 
29330 












32500 


1.241 


0.005 


0.109 


1.5 
2.0 


32050 
31580 






















35000 


1.539 


0.004 


0.103 


1.5 
2.0 


34540 
34060 























For A(refr) and T^ff the temperature Teg (y) is simply the elfective temperature written in the first column of the table. 
Parameters and ^j,^ are values of 5uv and 6ir calculated with T^ffiy) and fluxes FAlTcsiy), logg], which do not undergo the transforma- 
tion given by Eq.[6] To calculate them we have used - 1380 A and /i^ = 1 1084 A. Note that = Teff (y =1). 



2) Interpolate the model fluxes in the far-UV and IR for the 
adopted (J^s, log^) parameters; 

3) Use a least square procedure to search for an enhance- 
ment parameter j of the line blocking in the 1400 <A< 2150 
A wavelength interval, and correct the failure of the wind-free 
model fluxes used to fit the far-UV spectral region. The validity 
of y is then extended to the entire 0</i<2150A spectral re- 
gion. The empirical method used to calculate y and to modify 
the line-blocking is explained in ^3.2.1t 

4) Define a new effective temperature, reff(7'), to account 
for the lack of energy in the far-UV, produced by the increased 
absorption induced by y and for the redistribution of this en- 
ergy in the longer wavelengths. The calculation of Titsij) is 
explained in 33.2.2t 

5) Interpolate the model fluxes in the far-UV and IR for the 
reff(r) and logg = log§[reff(r),je] parameters; 

6) Calculate the angular diameter 6^ using Eq. (|2| and the 
IR fluxes dependent on the \T^s(j),\Qgg(J^'s(y),p)\ pair of 
fundamental parameters. With the fluxes interpolated in step 
5), calculate also the bolometric corrections 5uv and as 
indicated in Eq. ([3]). The bolometric correction (Juv is calcu- 
lated using the y-modified far-UV fluxes Fx defined in ^3.2. 11 
by Eq. where T^s is replaced by TisijY 



JUV 



7) Use i5uv and (5ir, and 6^ derived in step 6) to calculate 
the bolometric flux / with relation (O and obtain from Eq. ([1) 



,2 JO 



fobs 



(4) 



a new estimate of r{,.; 

eit 

8) Continue the iteration in step 3) by searching for a new 
enhancement parameter y using the far-UV fluxes Fx[Tf.ff(y)] 
interpolated in step 5). 



3.2.1 . Modification of theoretical fluxes 



Since we are not interested here in reproducing detailed en- 
ergy distributions to fit the observed ones, but rather in esti- 
mating the integrated amount of the emitted energy over the 
whole spectrum, we use an empirical method to minimize the 
disagreement between the predicted and the observed far-UV 
energy distributions. However, the method is not able to sim- 
ulate the consequences due to the presence of winds which 
become conspicuous in the extreme- UV (ISmifliet al.ll2002h . 
Nevertheless, this fact cannot significantly change the estimate 
of r^jj,, as we shall see in ^5.11 

Thus, we 
have proceeded in a sim ilar way as previously attempted by 
Remie & Lamersl(ll98lb and lZorec & Mercado-Ibanez.(.1987.) . 
i.e. we have modified the blocking degree of spectral lines in 
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the A ^2150 A region by writing the radiative flux emitted by 
a star at a given A as: 



(5) 



where Fc is the continuum flux for a given set of pa- 
rameters (Teff, logg), and is the Hne blocking fac- 
tor The line blocking factor bA easily can be calculated 
the Fa and F^ fluxes listed by CastelU & Kurucz 



using tne F^ and F,. fluxes listed by 
(l2003b (i http://kurucz.harvard.edu/grids.html| i. The enhanced 
line blocking factor was then calculated by multiplying b^ by 
a parameter y, which we assumed constant over all the wave- 
length interval < /i < 2150 A. The modified theoretical flux at 
a given A is then: 



FAiT^ff, log g) = FAiTes,logg){^YZy^) 



(6) 



For each star and at each iteration step of its T^^, we looked 
for the value of y that produced the best possible fit between 
the observed and theoretical fluxes, in the 1400 ^ A ^ 2150 
A wavelength interval. The model fluxes are calculated for a 
new effective temperature Tef{(y), which is different from T^^, 
and for a log g, which change as the iteration of the effective 
temperature continues. 



3.2.2. Effective temperature TeflfCy) of models 

It can readily be understood that the estimate of 5uv with the 
modified fluxes given by Eq. where in most cases y > 1, 
leads to an effective temperature lower than the nominal value, 
Teff, of the model used. Thus, when using Eq. ^ in the wave- 
length range 0</l<2150A, wecan recover the bolometric flux 
initially represented by the model characterized by {T^s, log g), 
by means of a new effective temperature, named hereafter 
Teff (y), which is larger than Teff if y > 1 . Mathematically, this 
is due to the lack of energy produced by an increased absorp- 
tion in the far-UV. Physically, this accounts for a redistribution 
of the excess of absorbed energy in the far- and extreme-UV, 
towards the near-UV, visible and IR spectral regions produced 
by the back- warming induced by the enhanced line blocking. 
To find the relation between Teif and Teif (y), we write the same 
bolometric flux, Fho\, in terms of both temperatures. On the one 
hand, it can be calculated using model fluxes dependent on Teff 
as: 



noi= FA{T,ff)dA+ FA{T,fi)dA, 

Jo JAa 



(7) 



and, on the other hand, the same Fboi can be obtained with mod- 
els that are a function of Te(f(y), as follows: 



noi= ^'.i[reif(y)]d^ + 
Jo 



/-•CO 



[reff(y)]d^, 



where /^,i[reff(y)] is given by the relation (|6|l. Since 



^bol - — -feif 

n 



(8) 



(9) 



and 

O-R 



-Tt^iy)^ FA[T,n{y)]dA+ F,[reff(y)]d^, (10) 
^ Jo Ja„ 

by subtracting (O from ^ and this result from (fTOb . we obtain 
the sought relation: 

T^siyf = + {—) \FA[T,s(y)] - FA[T,«{y)]]dA , (11) 
o"« Jo 

where we have not made explicit the dependence on log g, but 
have written that in the ^ /i ^ 2150 A wavelength interval, 
the unmodified fluxes Fa and modified fluxes Fa, are both cal- 
culated for the effective temperature reff(y). It is obvious that 
to have the coiTesponding reff(y) at each iteration step of the 
stellar effective temperature we put Teff = T^^^ and iterate the 
relation ( fTTT i. Then, replacing Eq. ^ into ( fTTT ). we derive: 

rfff = reff(y){l - (y - l)A[reff(y), logg]\" , (12) 
where the function A[retf(y), logg} is given by: 

A[reff(y),logg] = , '^ dA , (13) 

which can be calculated as a function of (T^sAogg) and used 
to derive the required value of A[reif(y)] by interpolation. The 
function A[reff(y), logg] is given in Table [1] where we also 
give the values of T^g, derived from ( fT2l i. In this table Tf^g{y) 
appears as the entering temperature. Actually, in the calcu- 
lation of the stellar effective temperatures, we enter relation 
(fT2l i with r^jj and deduce the corresponding value of Tf^^iy) 
at each iteration step. At each iteration step of T obtained 
by means of Eqs. ([T]) and and (Juv calculated with Q, we 
changed the gravity paramet er accordingly using the tables of 
log g - log gifi, Teff) given by Castelli & Kurucz ( 20061) . For in- 
stance, in Fig.[T]we show the results obtained at the final step of 
the analysis carried out for HD 41 1 17. For this particular star 
we have obtained Teff = 20740 K, refF(y) = 21 100 K, y= 1.53 
and logg = 2.5. All fluxes in FiglU observed and modelled, are 
normalized to a given flux in the visible wavelengths where the 
angular diameter is calculated. However, in the plot the loga- 
rithm of the fluxes is shifted by a constant value. We can also 
1.0 the model fluxes in the far-UV 



see that for T^^ and y 



are higher than the observed ones. In spite of the roughness of 
our approach, the fluxes from models computed with reff(y) 
and y - 1.53 improve the fit of the observed energy distribu- 
tion. It is also seen that the y-modified flux is slightly higher 
in A > 2200 A than the fluxes calculated with T^^, which leads 
to a smaller Balmer discontinuity as desired, i.e. D^bs - 0.050 
dex, D[reff(y)] = 0.065 dex, while it is DiT^^) = 0.074 dex. 

It is important to note that y > 1 produces a lowering of 
fluxes in the far- and extreme-UV, while the effect carried by 
the stellar wind on the photosphere , neglec ted here, increases 
the emitted fluxes (see ISmith et al.l (120021) V As we shall see 
in ^5.11 this increase happens at global flux levels that may 
have an effect on the estimate of 5uv at effective temperatures 
higher than 20 000 K, as shown in ^5.11 In general, the values 
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of 5uv obtained with wind-free models can differ significantly 
from those derived with non-LTE BW models if y > 1.5. The 
variation of 5uv with y can be obtained by noting that we can 
write: 



5uv(r) 



(14) 



where is given by Eq. We introduce the bolometric cor 
rections (5Jjy and (5Jjj defined as: 



Sum - 



UV ~ ^Ai, ' ^IR ~ ^Ai, 



(15) 



which are written in terms of fluxes calculated for reff(y) that 
do not undergo the transformation used for in the far- and 
extreme-UV. Subtracting 6^^^ from ( fT4b and taking into account 
the definition of Tss(y) written as: 



o-R, 



nAt 



n^ij)^ F,[reff(y)]d^x(l + 



by making use of (fTSl l. we readily obtain: 



(16) 



(^uv(r) = 5uv-(r- 1)(1+<5uv+^ir) X MT,B(7),logg] , (17) 

where (5Jjy, 6'^^ and the function A are given in Table [1] 

In the visible wavelengths used to calculate the stellar an- 
gular diameter, the blocking factor is ~ and the enhance- 
ment parameter is reduced to y = 1 . However, the model fluxes 
employed to obtain 0,i in Eq. (|2]i are for the effective tempera- 
ture Teff (y), i.e. F^ [7"eif (y)] ■ In general, when y 1 the effective 
temperature issued from ([T]i is higher than for y = 1 . Figure 



|2] shows Ar^jj - T^giy 1) - T^^^iy =1) for our sample of 
supergiants. Thus, even though y can sometimes be as high as 
y ~ 2, we see that the differences in the final value of T^^ are less 
than 300 K, which is much lower than the average expected er- 
ror affecting the determination of rfn, (see Finally, we note 
that for the stars in common with Remie & Lamers ( 198 1 ), our 
y values are systematically smaller by Ay ~ 0.5 than theirs, 
probably because in the models used here the line blocking ef- 
fect is more realistic than in the Kurucz (1979) models. 

4. The program stars and the observed quantities 

4.1. BCD parameters 

The program B-type stars, dwarfs to supergiants, are simply 
those for which both the BCD {Ai,D) parameters and calibrated 
fluxes from the far-UV to the IR, at least up to lfj.m were avail- 
able. In this work we excluded stars with the Be phenomenon, 
but included some AO to A2 type stars in the cold extreme of 
the hot fold of the BCD {Ai,D) diagram. The list of the program 
stars and their {Ai,D) are given in Table |2l 

Our sample contains 217 stars with MK luminosity classes 
from V to la. Observations in the BCD system were carried 
out over different periods. Most of them were observed with 



rf 



5 .• »vc 



• # 
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Fig. 2. Difference Arf^ = ^^.(y ?i 1) - ^^.(y = 1) against y. 



the Chalonge spectrograph te aillet et al.l 1973 1) attached to sev- 
eral telescopes at the Haute Provence Observatory from 1977 
to 1987 and at the ESO (La Sifla) from 1978 to 1988. More 
recently, on January 2006, low resolution spectra in the optical 
wavelength range AA 3400-5400 A were obtained at CASLEO 
(Argentina), with the Bofler & Chivens spectrograph. These 
spectra were wavelength and flux calibrated. Since the param- 
eter Al depends on the spectral resolution, we selected a reso- 
lution of 7 A at /I 3760 A, which is similar to the resolution re- 
quired in the original BCD system. From the spectra obtained at 
CASLEO, the BCD (Ai and D) parameters were directly mea- 
sured on the spectrograms. We have determined the parame- 
ter D = l og(F+37oo/F_37oo) followi ng the original BCD pre- 
scriptions (iChalonge & Divanlll952l) . We also controlled that 
^^-3700 corresponds to the flux level where the higher members 
of Balmer lines merge. Most stars have been observed many 
times, so that for each star the adopted (Ai,D) set comes from 
2 to about 50 determinations. The uncertainties that character- 
ize these quantities are then cr^, ^ 0.015 dex and cr^j ^ 2A, 
respectively. 

4.2. ASED data 



The ASED data were collected in the CDS astronomical 
database. The far-UV fluxes used in this work were obtained 
with the lUE spectra in low resolution mode and the 59 narrow- 
band fluxes, between 1380 and 2500 A m easured during the 
S2/68 experiment from t he TDl sateUite I jamar et af]|l976l: 
Macau-Hercot et al.lll978D . We have compared the TDl fluxes 
with the low resolution lUE spectra calibrated in absolute 
fluxes, but we have not detected any systematic deviation nei- 
ther in the far-UV ASED nor in the values of / obtained. The 
TDl fluxes give a homogeneous far-UV flux data set, for they 
are laboratory based calibrations. 

In the visible a nd near-IR spect ral region, fluxes are from 
Breger's catalogue ( Brepeij[T976alb) and the 13-color photom- 
etry calibrated in absolute fluxes dJohnson & Mitchelllll975h 
(hereinafter JM). The norm alized 13-color fluxes of the stars in 
common were compared to Breger ( 1976allbl) spectrophotome- 
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try and to the monochromatic fluxes observed by lTiig (Il98d) . 
The comparison of JM's absolute fluxes with those of known 
flux standards a Lyr (HD 17216 7), 109 Vir (HP 130109 ) and 
T/UMa (HD 120315) given in the lHaves & Lafliaml ( ll975b svs- 
tem revealed no noticeable differences in the /1/10.58-0.80 yum 
spectral region, which was chosen in this work to obtain the 
stellar angular diameter 6^ . Nevertheless, there are small dif- 
ferences near the BD that lie within the uncertainties of other 
calibrations, but they do not affect the estimate of /. We also 
noticed that JM's fluxes give consis tent c ontinuations to the 
lUE and ANS fluxes JWessehus et alJll982l) . 

4.3. Correction for interstellar extinction 



The adopted ISM extin ction law is from^ Cardelli et al" ( 1989h 
and lO'Donnelll ( Il993l) . The ratio = Ay/EiB - V) 
was adopted in this work accor di ng to the galactic re 



. K F 

gion, as specified in IGulati et alJ (119871) and IGulati et al 



(119891) . The adopted color excess E(B — V) is the average 



of several more or less independent determinations based on 
the following methods: a) UB V pho t ometry with the stan- 
dard intrinsic colors given by iLand (Il992h : b) BCD sys- 
tem with the intrinsic gradients O"^. tak en from the calibra- 
tion done by Chalonge & DivarJ ( 1973 ) from where we ob- 
tain that E(B - y) = 0.55 [ (brb - 0°J; c) depths of the 
2200 A ISM absorption ban d dBeeckmans & Hubert-Delplace 
1980t IZorec & Briot' 1985'); d) profile parameters of the 
2200 A band ( Guertler et a l. 1982; Friedemann et al. 1983; 
Friedemann & RoedeJl987h cahbrated in E(B-V) (iMoujtahid 



1993h : e) diagrams of E(B - V) vs. distance obtained with 'nor- 
mal' stars suiTounding the program stars within less than 2". 
The UBV photometry used is from the CDS compilation and 
the distances are either from the hipparcos satellite or spec- 
troscopy. 

4.4. Results 

The list of the program stars is given in Table |2] 1) HD num- 
ber of the star; 2) the MK SpT/LC determined with the BCD 
system; 3) the parameter Ai in Ai - 3700 A; 4) the Balmer 
discontinuity D in dex; 5) the adopted average color excess 
E(B - y), in mag; 6) the effective temperature T^^ with its esti- 
mated uncertainty, given in K; 7) the angular diameter with its 
estimated uncertainty, given in mas (milliarcseconds); 8) the 
parameter y used to fit the observed far-UV ASED (see ^3.2l l; 
9) the Tsfi(Ai,D) read on the new cahbration curves (see Fig. 

my 

5. Comments on the uncertainties affecting tlie T^^ 
and 9f determinations 

5.1. Systematic deviations 

Hot supergiants have massive winds whose optical depths can 
be high enough to heat somewhat the photosphere back to the 
continuum formation region. Thus, both spectral lines and the 
emitt ed continuum energy distr i bution can in princ i ple be mod- 
ified dAbbott & Hummerlll985l;[Gabler et al.lll989HSmith et al 



2002 : Morisset et al.ll2004 ). Since the models used in this work 
are wind-free, doubts can be raised about whether the use of 
these models introduces systematic effects on the estimate of 
T^g. As is seen in ^ this question may be of particular in- 
terest at Teff - 25 000 K. We then compared the absolute en- 
ergy distribution produ ced by a non-LTE BW model atmo- 
sphere for supergiants (ISmith et al.ll2002h with that predicted 
by an LTE wind-free model atmosphere for Teg = 25 000 K 
and logg = 2.95, similar to those used in the present work 
dCastelU & Kurucj |2003 ) . This comparison is shown in Fig. [3] 
where we see [Fig. [3ja)] that the most significant differences 
appear only in the extreme-UV energy distribution, at A ^ 500 
A. In this spectral region the fluxes are between two and three 
orders of magnitude lower than in slightly longer wavelengths, 
where the most significant contribution to the value of the the 
filling factor 6\jv in relation (|3]l arises. Although the differences 
seen in A ^ 500 A can certainly be important for the excita- 
tion/ionization of the stellar environment, as we shall see be- 
low, they seem to have a limited incidence on the estimate of 
the stellar bolometric flux / on which the estimate of T^^ relies. 

In Fig. [3b) we also see that the fluxes produced by both 
types of models are similar in the visible wavelengths where the 
angular diameter 0^ is calculated. This ensures that 6-^ does not 
suffer from the use of plane-parallel model atmospheres either. 
The same conclusion is reached from Fig.|5]by comparing the 
obtained 0^ values with the measured angular diameters. 

Let us now explore what systematic deviation can be ex- 
pected in our T^^ estimates due to the use of fluxes predicted 
by wind-free models of stellar atmospheres. In what follows, 
we use the notation T^^^ and 6'^^ to designate the parameters 
that should be derived with non-LTE BW model atmospheres. 
In ^6. H and Fig.|5]it is demonstrated that the angular diameters 
6^ obtained agree well with the observed ones. This means that 
they cannot introduce a systematic deviation in the derivation 
of T^g with Eq. ([T]i. Then, the ratio between T^^"^ and T^^^ is 
given by: 



eif 



eff 



1 + AW + AW 
^ ^ "UV ^ "iR 

1 + ^uv(r) + <5iR(r) 



(18) 



where (Juv(t) and 5iR(y) are the bolometric corrections cal- 
culated with wind-free models for different values of the en- 
hancement parameter y. The value of 5uv(t) is obtained from 
Eq. ( fTTb . while it is reasonable to assume that ^J'r -i5ir(7) = i5jj^. 
Concerning the non-LTE BW model (5[^y bolometric correc- 
tion, we ask what extreme-UV flux excess carried by the wind- 
related effects must exist to explain the underestimations sug- 
gested by the comparison shown in Fig.|2]that possibly affect 
our T^r,. determinations for supergiants. To this end we simply 



write: 



6vv(w) -wx 5jjv 



(19) 



where the factor w mimics the extreme-UV flux excess due to 
wind-related effects. Using for and A(reif , logg) the val- 

ues given in Table [1] we readily obtain the estimates of w for a 
series of imposed underestimations AT^^ which are displayed 
in Tabled 
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Smith et al. (2002) 
Castelli & Kurucz (2003) 



(b) 



Smith et al. (2002) 
Castelli & Kurucz (2003) 



A(A) 
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7000 
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Fig. 3. Comparison of theoretical ASEDs from non-LTE blanketed model atmospheres with stellar winds for T^g - 25 000 K and 
logg - 2.95 calculated by Smith et al. (2002) (doted line), with LTE wind-free models obtained by Castelli & Kurucz (2003) 
(full Une). 



Table 3. Extreme-UV flux excess factors w as a function of the 
efifective temperature and y, needed to explain the AT^^ under- 
estimations affecting the T^^. values for supergiants. 



K K 




y=1.0 7=1.5 7 = 2.0 


10000 -500 
-1000 
-1500 
-2000 


10.73 9.96 9.19 
21.94 21.02 20.09 
34.80 33.69 32.59 
49.45 48.14 46.83 


15000 -500 
-1000 
-1500 
-2000 


1.78 1.36 0.92 
2.68 2.18 1.68 
3.64 3.08 2.52 
4.70 4.06 3.43 


20000 -500 
-1000 
-1500 
-2000 


1.32 1.05 0.79 
1.67 1.37 1.08 
2.04 1.71 1.39 
2.43 2.08 1.72 


25000 -500 
-1000 
-1500 
-2000 


1.20 1.06 0.92 
1.41 1.27 1.11 
1.64 1.48 1.32 
1.89 1.71 1.53 


30000 -500 
-1000 
-1500 
-2000 


1.14 1.04 0.95 
1.28 1.18 1.08 
1.43 1.32 1.22 
1.59 1.47 1.36 



The estimates are done for log g = 3.0 



From the non-LTE BW models published by ISmith et al 



(120021) we obtain, 



w(25 000, 2.95)= 1.05 
^(30 200, 3.14)^ 1.08 



(20) 



for Teff = 25 000 K, log^ = 2.95 and T^ff = 32200 K, logg = 
3.14,respectively, where Teff = 25 000 K corresponds to the 



effective temperature where ther e are strong deviati ons in the 
diagram of Fig. |7] Unfortunately ISmith et al.l (120021) have not 
made available fluxes for T^ff < 25 000 K. Then, by extrapola- 
tion and approximate c alculation with our c odes for extended 
spherical atmospheres (Cruza do et al.ll2007l) . we get the value 
w(20000, 3.0) - 1.04. The values of w derived here are close 
to those in Table |3] for y = 1.5. Since for the superginats that 
are in common with those in the t'^ category (see ^6.2l l. we 

have y < 1.5, it means that we may expect our T^^ to be sys- 
tematically smaller by less than 470 K, 480 K and 640 K, for 
objects whose temperatures are 20 000 K, 25 000 K and 30 000 
K, respectively. In addition, for later B sub-spectral types, the 
underestimates can be even smaller, as it is difficult to believe 
that the w parameters are larger than those quoted above when 
TefF < 20000 K. For 7^2 (see Table |3]l differences easily can 
approach 1 000 K. Nevertheless, there is only one common late 
type supergiant with y > 1 .5 (HD 202850, y - 1 .78), for which 
curiously the difference between our estimate and that in the 
^eff gio'^P is +170 K. We note that to obtain the imposed un- 
derestimation Ar^jj. when w < 1 for a given factor y, the non- 
LTE BW models should predict lower extreme-UV fluxes than 
wind-free models do. 

The model atmospheres used in this work to obtain T^^ are 
interpolated here for log g parameters that were estimated us- 
ing the uvby-y6 photometry. As the index /3 is calibrated mainly 
for dwarf and giant stars, to see whether its extrapolation in- 
duces significant systematic effect for supergiants, we compare 
in Fig. |4]the log g(J3) adopted in the present work, with the 
log g(lines) parameters derived in the literature by detailed fit- 
ting of spectral lin es with model atmospheres, mainly hydrogen 



Hy and U6 fines (McErlean et al. 1999. 'Repolust et all 



Martins et al.ll 2005c iCrowflier et al.|2006;.Benagfia et al 



Markova & Pu^ l2008l; ISearle et al.ll20oi) . From this compari 



2004 



2007 



son we can conclude that our estimates of logg do not deviate 
either strongly or in a systematic way from those based on spec- 
tral lines. Although in the cited works the authors claim that 
their logg(lines) are determined with errors ranging from 0.10 
to 0.15 dex, they differ among them by 0.10 to 0.29 dex. The 
comparison in Fig. |4]of our log g(J5) estimates with log g(lines) 



10 



J. Zorec et al.: Fundamental parameters of B Supergiants from the BCD System 



shows that most of the points deviate randomly from the first 
diagonal within 0.25 dex. This implies that we do not expect 
that the log g(J3) values used in the present work will induce 
systematic errors in the T^^ and 6-^ parameters. 

In conclusion, we do not see how we can attribute to our 
BFM method systematic effective temperature deviations at- 
taining 2 000 K, or more than 5 000 K, as observed in Fig.|7] 



5.2. Random errors 

The T^g values issued from relation ([T]l and the 9-^ parameters 
derived with Eq. (|2]i have the following sources of error: a) the 
ISM color excess E(B - V); b) the line blocking enhancement 
parameter y; c) the log g parameter on which depend the model 
fluxes used to estimate the angular diameter and the filling fac- 
tor due to the unobserved spectral region; d) the filling factor 6 
used to calculate the bolometric flux in relation (O. 

In Appendix ^we explain in detail how we have calcu- 
lated the uncertainties of T^^j. and 6-^ by taking into account the 
combined effect of all error sources mentioned above. Since 
the most important uncertainty on T^^ and 6-1 is produced by 
the error on the ISM color excess estimate E(B - V), in this 
section we consider only the effect caused by this parameter as 
if it were the sole error source. Let us recall that the E{B - V) 
affects the bolometric flux estimate through /obs in (O of ^3.1! 
and the monochromatic fluxes on which the calculation of the 
stellar angular diameter depends. However, from ([T) and (|2| we 
find that: 



Table 4. Uncertainties on tL and 6^ as a function of Teff and 



rf 

eft 



(///.) 



1/4 



(21) 



which implies that the effect on the estimate of T^^ by an error 
on E(B-V) is reduced by the error of the monochromatic fluxes 
entering the calculation of 6-^. In Table|4]are given the estimates 
of the errors produced on T^^ and O'^ by the uncertainty AE(B - 
V) in the ISM color excess, assuming that it is the only source 
of error. We note the asymmetric propagation of uncertainties 
in and 0^ due to those in E{B - V). Since for us AE{B -V) ^ 
0.03, from Table |4] we find that on average the corresponding 
errors are of 5% in T-C and 2%in0^ . 



6. Comparison with other G and T^ff determinations 

6.1. Apparent angular diameters 

In spite of the fact that for the wavelengths where 6^ is calcu- 
lated, all models predict roughly the same flux level, the angu- 
lar diameter is in principle the quantity among those estimated 
in this work that can be the most strongly model-dependent. To 
test our angular diameter determinations, we have compared 
them with those found in the literature. We used the data col- 



lected in iPasinetti Fracassini et al.l (1200 Ih . where we discarded 
the very old determinations and privileged those determined 
by interferometry. The comparison is shown in Fig. |5] where 
w e put on the ordinate axis the a verage values from the data 



' eft 



log g due to errors in the estimate of the ISM color excess E(B- 
V). 



K 



AE(B - V) mag 
-0.07 -0.03 0.03 



0.07 



AE(B - V) mag 
-0.07 -0.03 0.03 



0.07 













= 4.5 








10000 


-0.054 


-0.025 


0.027 
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Fig. 4. Comparison of log g(Ji) parameters used in this work 
and derived with the uvby-y6 photometry, with those for the 
common supergiant stars estimated by other authors through 
spectral line fitting with model atmospheres. 



m 



Pasinetti Fracassini et alJ (l200ll) and the respective 1 cr error- mas. The uncertainty that can aff^ect the values due to eiTors 



bars. In this figure we can see that there is no systematic de- 
viation of points from the first bisecting line and that they dis- 
tribute around it with a fairly uniform dispersion: erg = 0.045 



f f 

in the angular diameter estimates ranges then from AT^^/T^^ - 
1.5% to 9% as the angular diameter goes from = 3.0 mas to 
6 = 0.5 mas. 
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Fig. 5. Average apparent angular diameters {6} in mas of our 
program stars determined by other authors (ordinates) against 
6^ obtained in the present work (abscissa). 




Fig. 6. Effective temperatures of dwarfs and giants determined 
by other authors (ordinates) against the T^g estimates obtained 
in the present work (abscissa). 



6.2. Effective temperatures 

In order to test the accuracy of the effective temperatures ob- 
tained with the BFM method, we compared our values to those 
found in the literature. The T^ff values were gathered according 
to the method used to determine them: 
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0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 

Fig. 7. Effective temperatures of supergiant stars determined 
by other authors (ordinates) against the T^^^ estimates ob- 
tained in the present work (abscissa). The eiTor bars correspond 
to tem peratures inside the ellipse taken from ICrowther et al.l 
(l2006l) (verticaD and in the present work (horizontal). The 
square with a downward error bar indicates the systematic av- 
erage shift that the .McErlean et al.. (.1999.) data might have. 



1) The temperature is derived empirica lly using either th e 
reddening-free color index QUV, as in lGulati et al. I (ll989h . 
or t he Co and ind ices of the Stromgren uvby-yS photome- 
try (|Casteni[l99i|)- 

2) The temperature is determined from integration of fluxes 
over a large interval of wavelengths, a s was done by 



Code et al.1 ( fl976l) : IUnderhill et all (ll979l) : lMalagnini et al 
(Il986l) . 



3) The temperature is obt ained by means of the BFM by 



Remie & LamersI (119811) . 



4) The temperature is evaluated by comparing the observed 
visual or UV fluxes with predictions from line-blanketed 

LTE mode l atmospheres , 

as in Malagnini et al. ( 1983): Morossi & Malagninil (119851) : 
iMala gnini & MorossL(,1990) . 

5) The temperature is computed by means of: a) NLTE line- 
blanketed model analysis of ionization balances due to 
He :/ II, Si m/ iv and Si ii/ m on mo derate resolution spectra, 
as done bvlMcErlean et al] (1 19991) using the code TLUSTY 
(|Hubenylll988h . or b) LTE line-blanketed model fitting of 
the optical region and Hy profile using ATLAS9 model at- 
mospheres as in Adelman et al. (2002). 

6) The temperature is determined from synthetic fits to the 
optical spectral range of B supergiants employing uni- 
fied NLTE line- and non-LTE BW extended model atmo- 



sphere cod es, such as F ASTWIND i Puis et al.1 l2005h o r 
CMFGEN ( jHiUier et al.l 12003) , as in Isearle et all ( l2008l) : 



Markova & Puis. (2008) : Crowther et alj (l2006l) . 
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Fig. 8. Comparison between the original BCD Balmer discon- 
tinuities of the program stars (abscissa) and those determined 
with absolute flux-calibrated energy distributions (ordinates). 



For each program star, Table |5] lists our determinations of 
eff together with the corresponding uncertainties and, in the 
following columns, the 7"^^ (/ - 1,...,6) obtained by other au- 
thors with the methods described in items: 1) to 6). When more 
than one determination of T^^ exists for the same star, obtained 
using similar techniques, we adopted the mean value. 

The difficulties in determining T^s can be better shown by 
separating dwarfs and giants from supergiants. Figure |6] dis- 



plays the comparison of r*^L' (stars), T^t! (open circles), T 



(2) 



r(3) 

'- eff 



(open triangles), T*^? (diamonds), T^^' (open squares), and T^^^ 
(plus signs) with our T^^^ values (abscissa) for dwarfs and 
giants. From this figure, we can see that the differences in 
the Tetf estimates seem to be related to the absolute value 
of temperatures rather than to the adopted techniques. For 
Jeff ^ 20000 K, all methods produce temperatures within 
o-T,^ = 750 K, while for T^s > 20000 K, ctt,,, = 1 300 K. 
We notice, however, that the T^^ derived with the BFM agree 
best with t hose obtained using line-blanketed LTE model at- 
mospheres (Malagnini et all 1983 : Morossi & Malagninil 19851 



Malagnini & Morossilll99ol) 

The same type of comparison, but for supergiants, is shown 
in Fig. |2] which appears less ordered than the plot for dwarfs 
and giants. However, the jf^ values do deviate strongly and 
in a syste matic way. The T^^} va lues for supergiants are taken 



only from McErlean et al. ( 19991) . and they were obtained with 
non-LTE wind-free model atmospheres. The observed devia- 
tions form a kind of 'temperature-step' rising at Teif ~ 15 000 
K and Teff ~ 25 000 K, i.e. at te mperatures identified as specific 



to the bi-stability phen omenon (Lamersetal 



1999t ICrowther etal] 12006.: Markova & Pulsll2008h . This fact 



1995t Vinketal 



stellar winds which have not been taken into account in their 
models. Moreover, the noted deviations also reveal a depen- 
dence on the method used to det ermine th e effective temper- 
ature. The Teff values from McE rlean et al.l (1999) determined 
by He n line profile fits deviate on average from our tempera- 
tures by (Areff ) = 5 50Qtl 800 K at abscissa <reff ) = 24 40Q+S00 
K; the temperatures determined by the Si iii/Si iv ionization bal- 
ance deviate by (A^eff ) = 1 300±380 K at {T^s) = 20 400±1 900 
K; those determined from the Si ii/Si iii ioniation balance de- 
viate by (AFeff) = 2 200+ 1 400 K at {T^s} = 1 4600 ± 1 500 K. 
The Teff va lues fromlMcErlean et al.l (Il999h have random errors 



±1 000 K. iMcErlean et al.1 (Il999h pointed out that their effec- 
tive temperatures could in fact be 10% lower, which means that 
for the stars lying in the ellipse of Fig. |7] the temperatures are 



on average 3 000 K too high. Thus, excluding the T^^* values, 
the average dispersions with the remaining sources cited in Fig. 
Hare cr^^,^ = 1 250 K for Teff < 20 000 K, and ctt,,, = 2 800 K 
forTeff >20000K. 

Stars in the jf^ group are from several sources, but 
their temperatures were derived using the same models, al- 
though corresp o nding to different implementation generations. 



ICrowther et al.l (l2006t) note that the cuiTent non-LTE BW 
model atmospheres lead to lower effective temperatures by 
1 000 to 2 000 K, as compared to some e arlier determinatio ns in 
Xudritzki et al. (1999); Crowtheretall (2002); Repolust etaP 
(120041) . and that their latest effective temperature estimates 
have uncertainties of about 1 000 K. However, three of the stars 
lying in the ellipse of Fig. |7] HD 30614, HD 37128 and HD 
38771, were assigned temperatures that are from 2 300 K to 
3 700 K higher than those obtained in the present work. From 
the discussion in Sect. 15. H it appears that our method of deter- 
mining the effective temperature does not introduce systematic 
deviations larger than 500 K to 700 K as the effective tempera- 
ture goes from 20000 K to 30000 K. We suspect then that the 
spectral lines used to derive some Tf^ may still undergo a per- 
turbation that the existing non-LTE BW model atmospheres do 
not account for entirely. Among the phenomena that may ex- 
plain the pertu rbed lines are: non-isothermal m ulticomponent 
plasma effects JSpringmann & Pauldrachiri992 ), shocks and/or 
the presence of exo-photospheric density clumps. Regarding 
the latter, it is worth noting that the use of an average continu- 
ous stellar wind to mask a possible clu mpy environment, w here 



a mass-loaded wind must appear (Har tquist et al.lll986h . can 



lead to misleading conclusions. The radiation transfer effects 
expected from the resulting average opacity are quite d ifferent 
from those the clumpy environment is able to produce (iBoisse 
19901) . Furthermore, the actual opacity of the medium is higher 



may reveal that there can be effects related to the presence of 



than that expected from line diagnostic and the back warming 
on the photospheric layers could have different characteristics 
than those expected from winds with regular density and tem- 
perature structures. 

7. Calibration of the (/li,D) parameters into T^s 

7.1. The scale of the Balmer discontinuities 

In the original BCD system, the Balmer discontinuities were 
obtained by comparing newly observed stars with stars ob- 



J. Zorec et al.: Fundamental parameters of B Supergiants from the BCD System 



13 



2.5 



1.5 



1.0 



\ 

\ 

\ 








-JZ..J /A 


r\ 






<Ai>= 


=32.1 A 


'A 










1 Q>^o 

X. >^ 










^ \ 


\ 








V 

•\ 










\ 










• 




















(a) 











).0 0.1 0.2 0.3 0.4 0.5 

D (dex) 




0.2 0.25 0.3 

D (dex) 



4 0.15 0.2 0.25 0.3 0.35 0.4 0.45 

D (dex) 



Fig. 9. The iteration of = constant curves in the plane (Ai , D). a) T-'^^ - T^^{D) relations for stars in strips of total width h.A\ = 10 
A and average </li-3700>= 32.1 and 52.3 A. b) Bisecting ^AxiD) curve for 7^= 15 000 K (full line), between the Ai ^Ax(D\Tk) 
obtained from panel a) (dashed curve) and the similar one derived by a least squares fit of a polynomial passing through stars 
in the plane (A\,D) having effective temperatures Tk+ATk (dotted curve), c) Iterated correction of the Teg - Teff{Ai , D) - const. 
layout. {A\,D*) - coordinates of a test star • = abscissa of intersections between the A\ - constant line and the Tk - constant 
curves at iteration step "n"; □ - abscissa of intersections between the A\ - constant Une and the new - constant curves at 
iteration step "n+1"; ★ = program stars. 



served simultaneously for which D was known. The 'zero' 
of these D values was known within an uncertainty of some 
0.012 dex (L. Divan, private communication). Conversely, for 
our program stars, we have derived the BD using absolute cal- 
ibrated fluxes, and compared them with those determined in 
the original BCD system. The result is shown in Fig.|8] In this 
figure we can see that, on average, the deviation between both 
types of BD determinations amounts to the expected 0.012 dex 
in the region of early sub-spectral types, but that the difference 
is actually a function of D. The least-square fitted relation be- 
tween the 'old' and 'new' BDs valid for the 0.03 < D < 0.55 
dex interval is given by : 



: 0.032 + 0.817 X Doid + 0.524 x D^^ ■ 



0.775 xD^,^, (22) 



where the D are given in dex. Therefore, in this work we used 
the scale of BDs determined from absolute fluxes, which enable 
the measured Balmer discontinuities to be directly compared 
with models. All diagrams in the present work are also given in 
the flux-calibrated Z)- scale. 

7.2. The calibration 

A calibration of the {A\,D) parameters into effective temper- 
atures for dwarf to g iant B-type stars has been presented in 
Divan & Zored (Il982h . In the present work, we are mainly in- 
terested in extending the BCD calibration for B supergiants. 
To ensure a better consistency of the layout of iso-effective 
temperature curves, we need to calibrate the entire region in 
the bi-folded BCD {Ai,D) surface that corresponds to OB -type 
stars of all luminosity classes. So, the iso-effective temperature 
curves for stars earlier than A2 of all luminosity classes are 



obtained in two successive steps: 1) First, we obtain an approx- 
imate layout of curves T^s{A\,D) — Tu —const.; 2) The shape 
of the approximate curves — const, is then corrected by iter- 
ation. In what follows, a detailed explanation of the procedure 
used is presented. 

1) Approximate system of curves Teff (/li ,D) — Ty_ -const.: 

The system is established as the mean regression curve be- 
tween the variables x and y from a set of measured pairs of 
points {x,y). The regression x = jt:(y) does not necessarily repre- 
sent the same locus of points as the one calculated in the form 
y - y{x), nonetheless issued from the same (x, y) data set. The 
bisecting curve between the direct x = ^(y) and the inverse func- 
tion jc =y ' (jc) can then be used to represent the sought relation. 
Thus, in this work we calculate two series of functions to ob- 
tain the approximate average/bisecting system of A\ -A\{D\Tk) 
curvefl where Tk - const, with A: = 1, 2, ...A^k are constant val- 
ues of the effective temperature: 

/) A series of A^k regression polynomials A\ -A\{D) is cal- 
culated using stars whose effective temperatures are in inter- 
vals {Tk - ATk, Tk + ATk), where ATk/Tk = 0.05. For a given 
Tk the obtained regression polynomial is assigned to the aver- 
age Tk temperature calculated from those entering the interval 
(Tk - ATk, Tk + ATk); 

ii) With the T^g values of stars lying in successive strips 

of constant Ai parameters and total width AAi = 10 A in the 

f f 

{Ai,D) plane, we obtained a series of curves 7^^^ = ^^^^(DK/li)), 
where (/li) is the average of Ai parameters of the stars entering 
a given strip, as shown in Fig. [9^ for two (Ai) values. Then, 
for a given Tk temperature obtained in step /), we read the ab- 



' Here, the notation F = F(u\v) means that F is a function of u, so 
that v = v(F, u) = constant over the entire space of variables [F, u). 
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Fig. 10. (a): Spline-smoothed T^^{A\,U) - const curves for the values of effective temperature given in the box at right-top; 
(b) The curvilinear-quadrilateral BCD (X\ , D) spectral classification diagram, where the corresponding 2D MK spectral types 
(bottom shaded strip) and the luminosity classes (left shaded strip) are indicated. The cross-marks indicate the mid-point of each 
MK spectral type-luminosity class box where the Tf^f({A\,D) values given in Table 3 were obtained. 



scissae determined by all the - {Ai}(D) curves, as in the 
example shown in Fig. [9^ for the particular case Tk = 15 000 
K, to draw the corresponding {Ai} - {Ai}(D\Tii) (see Fig.|9j3, 
dashed curve). 

The final curve adopted as the first approach to Tii(A\ , D) - 
const., here called Ai = Ai{D\Tii), is the smoothed bisecting lo- 
cus of points (full line in Fig.|9}5) between the regression poly- 
nomial Al =Ai(D\Tk) obtained in /) (dotted line in Fig.|9]3) and 
the homologous one (Ai) - {Ai}(D\Tk) constructed in //). Using 
a second order interpolation, the set of Ai = Ai(D\Tk) curves 
with k - 1,2, ...Nk is then transformed into Ai = Ai(D\Tk) 
curves, where = 9 500, 10 000, 35 000 K are rounded in- 
teger values used in the following step. Hereafter the curves 
Ai(D\fk) are called Ai(D\Tk). 

2) Corrected Teff (/li , D) ^onst curves: 

The layout of Tk{A\,D) — const curves obtained in 1) in 
principle might be adopted as the sought calibration. We have 
noticed, however, that we can reduce even more the residu- 
als ej-^j = \T^s{A\,D*) - rfjfl, where T^s(A\,D*) represents the 
temperature read in the calibrated {A\,D) plane for the star 
whose parameters are {A*^,D* , T^^). To reduce e^^^ somewhat, 
we proceed as follows. We consider a star whose parameters 
are iA\,D*, T^^^). We call Dk-i, D^, D^+i and D^+i the abscis- 



sae of the intersections of the Ai - A*^ line with the curves of 
constant effective temperature Ai - Ai{D\Tk), ordered so that 
Tk-i > Tk > r^jj > Tk+i > Tk+2- We calculated the 'corrected' 
abscissae D^, D^^j of the curves Tk and Tk+i at Ai =A\ with: 

D' = D' + {2II^ll!I^)(Dk-i - D*) 
k li/r,-_,-i/rf/^ ^ ^23) 

For a given (/l*,D*), the abscissae Dk x and Dk+i are main- 
tained unchanged (see Fig. |9]:). This operation is done for all 
program stars, so that each of them contributes with its corre- 
sponding correcting displacement D' - D. In this way the po- 
sition of each curve Tk - const, depends also on the nearest 
ones. The smoothed curve passing through the new {Ai,D') 
points relative to the same Tk = const, is considered as the new 
Al = Ai(D\Tk) curve. Thus, at each iteration the whole sys- 
tem of Tk curves is somewhat modified. The iteration is contin- 
ued until two conditions are obeyed: 1) the correlation of the 
read effective temperatures of the program stars in the newly 
obtained layout of Teff - constant curves as a function of the 
respective T^„ values does not show any systematic deviation 
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from the first diagonal; 2) the dispersion of points around this 
line attains its lowest possible value. 

The Tf.ff{Ai,D) diagram thus obtained is shown in Fig.fTOk. 
The upper limiting curve in the diagram of Fig. [TOk corre- 
sponds to the lowest Ai parameters ever observed in the BCD 
system for OB supergiants. Whilst Ai values higher than the 
bottom limiting curve of this diagram have never been mea- 
sured for normal dwarf OB stars, they are normally observed 
for O sub-dwarfs and/or white dwarfs. 

The method used to obtain the curves of constant Tgff{Ai, D) 
makes the position of each T^ff - const curve depend on the 
global Tsff(Ai,D) pattern; i.e. a change in one of these curves 
has an incidence on the shape of the surrounding curves. We 
also notice that the layout of the Teff - const curves in the sector 
of supergiants is strongly constrained by the way the curves are 
arranged in the (/li,Z)) region of dwarfs and giants which are 
traced by a large sample of stars, and by the natural limiting 
condition which imposes that D — > as /li becomes negative. 
Actually, before the limit D — > is attained, for supergiants with 
Jeff S 25 000 K, there is a luminosity class interval, roughly 
from /li -3700 ~ 55 A to -3700- 5 A, the extended low 
density atmosphere of these stars behaves like the circumstellar 
envelope of Be stars, producing D < 0, i.e. an emission-like 
Balmer discontinuity. The black body-like behavior D = is 
reached for the theoretical limit Ai -3700 — 47 A. 

To examine the correspondence between the T^ff values and 
the MK spectral classes, we show in Fig.fTOb the calibration of 
the BCD (Ai,D) surface into B-spectral subtypes and luminos- 
ity classes. We can see that the pattern of T^ff = const curves 
mirrors the empirical MK spectral type calibration. The delim- 
iting strips of constant sub-spectral types are determined with 
color gradient curves O^i = const. This system of curves resem- 
bles each other because the color temperature in the wavelength 
region over which (^rh is defined is close to the stellar effective 
temperature. In Table|6]we give the values of T^g for the points 
marked in Fig.fTOb with crosses. These are mid-points of more 
or less large curvilinear boxes corresponding to a given MK 
SpT/LC. We note that it may not be straightforward to com- 
pare the effective temperatures for the nominal MK SpT/LCs 
given in Table |6] with the average values determined in the lit- 
erature for the same MK spectral classes. This is because the 
baricenters of the (/li, Z)) parameters in these averages may not 
correspond to the mid-points of the curvilinear quadrilaterals 
given here (see crosses in Fig.fTOb). 

Finally, to verify the reliability of the new BCD Teft calibra- 
tion, in Fig. [TT] we compare the effective temperatures of the 
program stars read in the {Ai,D) calibration to those derived 
with bolometric fluxes. The comparison reveals a good corre- 
lation, where the dispersion of points around the first diagonal 
remains within a constant limit AFetf /Teft - ±0.17. Considering 
that the T^^ values represent the reference for our Tf.ff{Ai,D) 
calibration, the read Tgff values in the {Ai,D) plane can be con- 
sidered as affected by two kinds of uncertainties, one of them 
inherent to the measurement of the Balmer discontinuity, and 
the other related to the obtained layout of the T^^ff - constant 
curves. Since the deviation of points around the first diagonal 
in Fig. [TT]has a Gaussian distribution, which lies in the whole 




tJ/io'* k 

Fig. 11. Comparison between effective temperatures read in the 
{Ai,D) plane (ordinates) and those derived with bolometric 
fluxes of the program stars (abscissa). 

3cr/7'eff = +0.17 interval, the expected total average erro r af- 
fecting the read Tefi(Ai,D) quantities is thus (ISmartlll958l) : 



ja^.±0.05, 

Tefl V 7T Tef[ 

which is valid for the entire calibrated {Ai,D) plane. 



(24) 



8. Discussion 

In this work we have obtained a new retf(/li,Z)) calibration 
which is homogeneous over the entire BCD plane of (/ii,Z)) 
parameters coiTesponding to early-type stars. For dwarfs and 
giants, this calibration yields Tefj values which are consistent 
with those derived using photometric and spectroscopic tech- 
niques, and with model atmospheres. However, for supergiants 
our Teif estimates show rather h igh systematic discrep ancies 
with the Tefl values derived by iMcErlean et al. I (Il999h with 
wind-free models of stellar atmospheres (see Fig.|7|l- As seen in 
Fig. I2I the wind-free model-dependent T^ff values are system- 
atically higher than ours. Our T^s estimates are also lower than 
those obtained with non-LTE BW models around T^ff ~25 000 
K, but the discrepancies are smaller. From the discussion in 
^5.11 it appears that our estimates for supergiants could 
have possible systematic underestimates ranging from 500 K 
to about 700 K as the effective temperature goes from 20 000 K 
to 30 000 K. Although this can partially account for the devi- 
ation of points belonging to the group (crosses), it cannot 
explain the points above the the one-to-one line in Fig.|2]from 
other groups, because the systematic underestimations of our 
r^jj. values are lower than 1 000 K. We suggest that the T^fi val- 
ues determined through line profile fitting with synthetic spec- 
tra obtained with non-LTE BW model atmospheres, which in 
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Table 6. Effective temperatures at the mid-point of each curvilinear quadrilateral representing a MK spectral type-luminosity 
class in the BCD (Ai,D) plane for early type stars. 



LC : 




V 




IV 


111 


II 


lb 


la 


SpT 


Al 


D 




Al 


D 




Al 


D 




Al 


D 




Al 


D 




Al 


D 




06-7 


67.0 


0.062 


34640 


59.4 


0.060 


34800 


54.5 


0.055 




50.2 


0.048 


33970 


45.9 


040 


32590 


37.3 


0.014 


29670 


08-9 


66.1 


0.075 


32750 


58.4 


0.071 


33070 


53.1 


0.067 


32460 


47.9 


0.061 


30690 


43.1 


0.052 


28350 


34.7 


0.025 


26040 


BO 


65.1 


0.094 


30000 


56.9 


0.091 


29970 


51.2 


0.085 


29070 


45.0 


0.077 


26710 


39.9 


0.066 


24380 


31.9 


0.037 


23890 


Bl 


63.5 


0.121 


27080 


54.9 


0.116 


26960 


48.4 


0.110 


24880 


41.5 


0.099 


23370 


36.3 


0.085 


22360 


28.8 


0.052 


22020 


B2 


61.5 


0.160 


22620 


52.2 


0.156 


22380 


44.9 


0.147 


21740 


37.5 


0.127 


20570 


32.2 


0.107 


19790 


25.0 


0.071 


19310 


B3 


59.7 


0.203 


19270 


49.7 


0.201 


19160 


41.6 


0.188 


18900 


33.5 


0.160 


18010 


28.2 


0.133 


17520 


21.6 


0.091 


17140 


B4 


58.8 


0.235 


17220 


48.3 


0.235 


17130 


39.6 


0.218 


17310 


31.1 


0.184 


16500 


25.7 


0.152 


16180 


18.9 


0.108 


15660 


B5 


57.9 


0.269 


15380 


47.1 


0.271 


15570 


37.8 


0.249 


15890 


29.0 


0.211 


15080 


23.1 


0.174 


14770 


16.6 


0.126 


14460 


B6 


57.7 


0.302 


14100 


46.7 


0.309 


14200 


36.3 


0.286 


14530 


27.1 


0.241 


13790 


20.9 


0.197 


13510 


14.6 


0.143 


13560 


B7 


58.1 


0.333 


13000 


47.0 


0.345 


13000 


35.3 


0.322 


13460 


25.7 


0.273 


12600 


19.2 


0.221 


12540 


12.8 


0.159 


12800 


B8 


59.1 


0.363 


12190 


47.7 


0.380 


12120 


35.0 


0.363 


12380 


24.8 


0.306 


11790 


17.9 


0.246 


11930 


11.2 


0.176 


12280 


B9 


61.5 


0.406 


11340 


49.5 


0.427 


11240 


35.8 


0.423 


11240 


23.9 


0.362 


10800 


16.0 


0.291 


11070 


8.9 


0.205 


11720 


AO 


64.8 


0.448 


10470 


52.1 


0.476 


10310 


37.3 


0.477 


10350 


23.5 


0.431 


10140 


14.3 


0.352 


10310 


6.3 


0.243 


11010 


Al 


67.4 


0.476 


9860 


54.0 


0.508 


9730 


38.4 


0.511 


9820 


23.7 


0.478 


9680 


13.4 


0.401 


9690 


4.3 


0.277 


10440 



LC = Luminosity class; SpT = spectral type 
[A,] =A, - 3700 A; [D] = dex 



most cases are isothermal, could be partially responsible for 
the observed effective temperature discrepancies. It should not 
be neglected either that some discrepancies could arise due 
to problematic convergences of mod els in the diluted exter- 
nal at mospheric layers of s upergiants ( Crivellari & Simonneaul 
19911: jshnonneau & Crive lM l l993l: 



19941 [simonneau&Crivellanlll994 . 



Crivellari & Simonneau 



8.1. Effects due to metallicity 

The present BCD reft calibration was obtained using a sam- 
ple of stars near the Sun, which are characterized by an av- 
erage metallicity Z - 0.02. We might then expect system- 
atic differences between our Tgff(Ai,D) estimates and those for 
stars having lower metallicities, as in the Magellanic Clouds. 
Since stars of the same mass, but with lower Z, have smaller 
radii, we cannot ensure that a given T^^. can correspond to 
the same pair (Ai , D) of a Magellanic Cloud-star. We can then 
wonder whether the T^jj-based calibration of the (Ai , D) pa- 
rameters obtained here can be used for stars with Z = 0.004 
and lower. To assess possible systematic differences in the 
Tes estimates caused by the use of the present reff(/li,Z)) cal- 
ibration for stars with an initial metallicity other than Z = 
0.02, we calculated the BD using model atmospheres for Z = 
0.002, 0.0002, and Z = 0.2, us ing the grids of ATLAS 9 
model atmospheres calculated by ICastelh & Kuruc j ( 2003 ). 
The results are displayed in Table [T] We have not extrapo- 
lated the fluxes to lower \ogg values because of the diffi- 
culties of model convergence that cannot ensure reliable D 
values . In this table we see that for given values of \ogg 
and Tf.ff, the D parameter is larger the lower the metallicity. 
Although the differences in the values of D depend on Teff 
and \ogg, on average we have D(Z = 0.002)-D(Z = 0.02) = 
0.014 + 0.003, D(Z = 0.0002)-D(Z = 0.02) = 0.023 + 0.005; 



D(Z = 0.2)-D(Z = 0.02) = -0.019 + 0.008. These differences 
are slightly larger than those expected from the intrinsic un- 
certainties of the D values. The present BCD calibration of 
effective temperatures would then produce slightly underes- 
timated Tfiff values for (/li,Z)) pairs of stars in environments 
with Z < 0.02, or overestimated if Z > 0.02. In both cases the 
Tgff{Ai,D) of stars in environments with Z 0.02 determined 
with the present calibration can be easily corrected using Table 
|7] The present r^^j-based calibration of the (Ai,D) may be in- 
teresting also for early-type stars in the Magellanic Clouds, 
since there are no systematic measurements of far-UV fluxes 



that would enable us to use the BFM to obtain their T 



In Fig. [12] we show the relation between the effective tem- 
perature and the Balmer discontinuity. Figure fT2k refers to 
the relations issued from the calibration obtained in this work, 
where the values are for the mid-points of the spectral type- 
luminosity class curvilinear quadrilaterals in the BCD plane 
of (Ai,D) parameters (c.f. Table |6l). In this figure we indicate 
the lines of constant mean MK spectral type of stars in dif- 
ferent average luminosity classes near the Sun. In Fig. fT2b 
we show the (Teff , D) curves calculated with wind-free mod- 
els dCastelli & Kurucal2003b for several \ogg values and the 
standard metallicity Z - 0.02. The gray-shaded region corre- 
sponds to (Teff , D) curves for Z - 0.0002 and for the same set 
of \ogg parameters as for Z = 0.02. For comparison sake, the 
hatched zone demarcates the region occupied by the empirical 
curves displayed in Fig. [T2h . This shows that while wind-free 
model atmospheres can account reasonably for the (Teff, D) re- 
lations from dwarfs to giants, they fail to do so for supergiants, 
and that the problem is particularly acute for supergiants with 
Teff > 20000 K. It is expected that non-LTE B W models of hot 
supergiants can explain more easily the observed (Teft, D) rela- 
tions, since such atmospheres are heated by the backward wind 
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0.0 0.1 0.2 0.3 0.4 0.5 0.0 0.1 0.2 0.3 0.4 0.5 

D (dex) D (dex) 



Fig. 12. (a): Effective temperature Tsft(Ai,D) against the Balmer discontinuity for different average MK luminosity classes and 
spectral types. The letters OBA head the lines of constant spectral type; (b) Teif = T^siD, \ogg) relations from wind-free model 
atmospheres with metallicity Z = 0.02. The gray region corresponds to Teff - Tf.f[(D, \ogg) relations for metallicity Z - 0.0002 
and the same log g parameters. The shaded zone corresponds to the empirical T^siAi , D) curves from panel (a) of this figure. 



radiation, and the corresponding D values should then become 
smaller. 



8.2. Effects related to the rotation 

The new calibration of the (Ai , D) plane as a function of T^s 
was done assuming that all stars can be c haracterized by two 
parameters: T^s and log g. It was shown by Fremat e t al.l (12005) 
that also the apparent visual spectral region emitted by the 
non-uniform and geometrically deformed surface of fast rotat- 
ing early type stars can be well represented with parent non- 
rotating Tetf and \ogg parametric counterparts. The present 
caUbration can then provide a first step for the interpreta- 
tion of spectra emitted by rotating earl y-type objects, as al - 



ready shown inlLevenhagen et al 
Zorec et al ](l2005h: Vinicius et al 



(2003):lNeiner et al 



(2006):lFrematet al 



(2003); 



( 20061) : 



Martavan et all (l2006il2007h : lFloquet et all (12000112003 ^ 

Rotation can also induce internal mixing of chemical ele- 
ments. The mixing and transport of chemical elements through- 
out the star, produced in the stellar core, can be stronger with 
higher stel lar rotation and it can be enhanced if the metallic- 
ity is low dMevnet & Maedei I2OOOI l2002h . In particular, this 
may concern helium, whose abundance could then be larger in 
the atmospheres of evolved stars. Since this excess of He con- 
tributes to the absorption in the stellar surface layers, the value 
of the Balm e r disc ontinuity can be affected. It was shown by 
Cidale et al. I (E007 ) that in hot dwarf stars the larger the He 
abundance the smaller the value of D. 

In order to see in more detail the influence of the He/H 
abundance ratio on the emitted visual energy distribution, a 
grid of synthetic spectra for a range of effective temperatures 



was computed in n on-LTE using the TLU STY and SYNSPEC 
computing codes ("Hubenv & Lanz and the references 

therein, assuming model atmospheres with He/H ratios of 0.1, 
0.2, 0.5 and 1.0, and Z - 0.02. The atomic models we used 
are basically those provided on the TLUSTY website for H i 
(9 levels). He i (20 individual levels) and He 11 (20 levels). In 
all cases, the micro-turbulence velocity was supposed to be 2 
km s We have reduced the synthetic spectra to the BCD res- 
olution and calculated the BDs following the procedure used 
for the empirical ones. In Table |9] we give the obtained differ- 
ences 6D - D-D(O.l) as a function of the model (Teff, logg, 
He/H) parameters. Here D(O.l) indicates that D is for the He/H 
= 0.1 abundance ratio. In this table we see that on average 
SD < and that \6D\ is larger the higher the He/H abun- 
dance ratio an d for 19000 < Tes < 2 3 000 K. In recen t studie s 
of supergian ts Repolust et al.1 ( 20051) : Crowther et al. ( 2006 ): 



Searle et al 



(120081) : iMarkova & Pulsl (I2008h have found that 
an abundance He/H — 0.2 fits their atmospheres. From Table 
|9]we see that a ratio He/H = 0.2 effects the value of D by less 
than the error of measurement of D. Nevertheless, these dif- 
ferences are systematic and can be important for the hottest 
supergiants. In table (Table |8]l we give the increments dT^a 
that can be produced on the Tetj values for log^ - 3.0 by 
5D{<d.2) = D(He/H = 0.2)-D(He/H = 0.1) with respect to the 
T^s{D) scale at He/H = 0.1 given in Table |7] 

We note that: 1) Only those stars that were fast rotators dur- 
ing their Main Sequence evolutionary phase are likely to show 
an atmospheric increase of the He/H ratio in their blue super- 
giant phase; 2) The energy distributions in the near-IR where 
the angular diameter 9^ is calculated, using models for He/H 
= 0.1 and He/H - 0.2, give the same flux levels. Moreover, in 
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Table 7. Balmer discontinuities D as a function of T^ff and log^ 
for different metallicities Z. 









lo 










4.5 


4.0 


3.5 


3.0 


2.5 


2.0 


Tea 






Z = 


0.02 






10000 


0.508 


0.520 


0.526 


0.521 


0.502 


0.452 


12500 


0.370 


0.361 


0.353 


0.342 


0.312 


0.281 


15000 


0.282 


0.270 


0.258 


0.237 


0.201 


0.168 


17500 


0.221 


0.209 


0.195 


0.169 


0.132 


0.097 


20000 


0.175 


0.165 


0.150 


0.121 


0.086 




22500 


0.139 


0.129 


0.113 


0.080 


0.048 




25000 


0.111 


0.101 


0.082 


0.052 


0.011 




27500 


0.091 


0.080 


0.061 


0.030 






30000 


0.072 


0.060 


0.035 


0.012 






32500 


0.055 


0.041 


0.027 








35000 


0.040 


0.027 


0.012 








Tea 






Z = 


0.002 






10000 


0.511 


0.523 


0.526 


0.518 


0.500 


0.457 


12500 


0.380 


0.376 


0.368 


0.357 


0.326 


0.295 


15000 


0.295 


0.288 


0.276 


0.257 


0.222 


0.190 


17500 


0.234 


0.227 


0.213 


0.188 


0.155 


0.122 


20000 


0.189 


0.181 


0.166 


0.139 


0.108 




22500 


0.153 


0.145 


0.129 


0.102 


0.074 




25000 


0.125 


0.115 


0.099 


0.073 


0.048 




27500 


0.101 


0.091 


0.074 


0.050 






30000 


0.082 


0.070 


0.052 


0.032 






32500 


0.065 


0.052 


0.034 








35000 


0.051 


0.036 


0.018 








Tea 


Z = 0.0002 


10000 


0.513 


0.524 


0.528 


0.521 


0.503 


0.462 


12500 


0.388 


0.385 


0.379 


0.366 


0.338 


0.308 


15000 


0.304 


0.298 


0.287 


0.266 


0.233 


0.201 


17500 


0.243 


0.236 


0.223 


0.199 


0.165 


0.133 


20000 


0.197 


0.189 


0.175 


0.150 


0.119 




22500 


0.160 


0.152 


0.138 


0.115 


0.087 




25000 


0.131 


0.122 


0.109 


0.088 


0.064 




27500 


0.106 


0.097 


0.084 


0.067 






30000 


0.086 


0.076 


0.064 


0.050 






32500 


0.068 


0.058 


0.047 








35000 


0.053 


0.042 


0.032 








T fr 
^ eii 






Z = 


0.2 






10000 


0.499 


0.515 


0.528 


0.525 


0.507 


0.458 


12500 


0.361 


0.357 


0.346 


0.334 


0.317 


0.299 


15000 


0.265 


0.258 


0.248 


0.222 


0.200 


0.181 


17500 


0.199 


0.190 


0.179 


0.151 


0.127 


0.104 


20000 


0.151 


0.142 


0.127 


0.103 


0.079 




22500 


0.115 


0.106 


0.089 


0.070 


0.048 




25000 


0.088 


0.079 


0.061 


0.046 


0.027 




27500 


0.067 


0.058 


0.041 


0.028 






30000 


0.050 


0.040 


0.027 


0.014 






32500 


0.037 


0.026 


0.017 








35000 


0.024 


0.014 


0.011 









Table 8. Increments dT^-s carried by 6D{Q.2). 









Tea (K) 






17000 


19000 


21000 


23000 


25000 


27000 30000 








STea (K) 






+140 


+ 170 


+270 


+590 


+600 


+610 +600 



changes were obtained in the value of the bolometric correc- 
tion 5uv between He/H = 0.1 and He/H - 0.2, so that the use 
of models for He/H - 0.2 in our BFM would not in principle 
change our T-^^ estimates; 3) Even though the calculations of 
5D given in Table|9]concern wind-free model atmospheres, the 
low dependence of deviations 6D with Teif and log^ ensures 
that the inclusion of the non-LTE BW will not change the re- 
sults noticeably. 

Table 9. Differences of Balmer discontinuities 5D = Z)(He/H)- 
D(O.l) at metallicity Z = 0.02, as a function of Teff, log^ and 
for different He/H abundance ratios. D(0. 1) is for the He/H= 0. 1 
ratio. 



He/H=0.2 0.5 1.0 



Tea 






5D (dex) 




12500 


3.0 


+0.011 


+0.007 


+0.000 




3.5 


-0.005 


-0.006 


-0.014 




4.0 


-0.001 


-0.004 


-0.008 


15000 


3.0 


+0.000 


-0.014 


-0.020 




3.5 


-0.005 


-0.014 


-0.029 




4.0 


-0.005 


-0.014 


-0.029 


17000 


3.0 


-0.003 


-0.019 


-0.025 




3.5 


-0.005 


-0.016 


-0.032 




4.0 


-0.006 


-0.016 


-0.035 


19000 


3.0 


-0.003 


-0.020 


-0.026 




3.5 


-0.004 


-0.016 


-0.032 




4.0 


-0.006 


-0.017 


-0.035 


21000 


3.0 


-0.006 


-0.019 


-0.025 




3.5 


-0.004 


-0.014 


-0.029 




4.0 


-0.006 


-0.016 


-0.033 


23000 


3.0 


-0.006 


-0.017 


-0.022 




3.5 


-0.003 


-0.013 


-0.026 




4.0 


-0.005 


-0.015 


-0.029 


25000 


3.0 


-0.006 


-0.014 


-0.019 




3.5 


-0.003 


-0.011 


-0.023 




4.0 


-0.004 


-0.013 


-0.025 


27000 


3.0 


-0.004 


-0.010 


-0.016 




3.5 


-0.002 


-0.009 


-0.019 




4.0 


-0.003 


-0.011 


-0.021 


30000 


3.0 


-0.003 


-0.005 


-0.011 




3.5 


-0.002 


-0.007 


-0.013 




4.0 


-0.002 


-0.008 


-0.014 



6D = Q for all He/H abundance ratios at Tea = 10000 K 



the very few cases with T^ff > 20 000 K where we have calcu- 
lated the energy distributions in the far- and extreme-UV, no 
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9. Conclusions 

In this work we have presented a new and homogeneous cal- 
ibration of the BCD plane (Ai,D) as a function of Tg^ for 
early-type stars of all luminosity classes, in particular super- 
giants, which complete a similar one made ea rlier only for 
dwarf to giant B-type stars dOivan & Zorecl[l982 ). The present 
calibration is based on effective temperatures calculated with 
the bolometric-flux method for all program stars, whose indi- 
vidual uncertainties are on average e/T^^^ - 0.05. The average 
error of the obtained Tes values on the {Ai,D) plane is the same 
in all early spectral types and luminosity classes, and they are 
of the same order as for the individual T^g values. The effec- 
tive temperatures of OB supergiants derived in this work agree 
within some 2000 K with other determinations found in the 
literature, except with those issued from wind-free non-LTE 
plane-parallel models of stellar atmospheres, which produce 
overestimates of up to more than 5000 K near Teff - 25000 
K. 

The Tgff{A] , D) calibration has the advantage of using mea- 
surable parameters of the continuum spectrum around the 
Balmer discontinuity, which are strongly sensitive to the ion- 
ization balance of the photosphere and to its gas pressure. The 
BCD parameters {Ai,D) can be easily measured in spectra 
of low resolution and, even though D is slightly affected by 
the ISM extinction, the correction can be precisely and eas- 
ily accounted for. This makes the {Ai,D) quantities useful for 
distant stars. Moreover, since their determination easily can 
be programmed for automatic measurements, they are conve- 
nient for studies of stars in clusters or other stellar systems ob- 
served with multi-object spectrographs and/or spectro-imaging 
devices. The BCD system can also be of interest in analyz- 
ing stars in far away galaxies that will be observed with the 
Extremely Large Telescopes in the near future. In this con- 
text, the present calibration of reff(/li, Z)), together with those 
in preparation for \og g{A\,D) and M\,o\{A\,D), can be help- 
ful in translating the data in terms of physical quantities such 
as Teff and log^. Furthermore, introducing small corrections, 
these calibrations can be applied to stars in environments with 
different initial metallicity. 

The contamination of atmospheres in supergiants by He 
due to rotational mixing does reduce the value of the Balmer 
discontinuity. However, this change can exceed the average 
measurement uncertainty of this parameter only when the 
abundance ratio becomes He/H>0.05, which is twice as large 
as the value assumed today in the model atmospheres of su- 
pergiants. Nevertheless, differences of up to 600 K can be pro- 
duced by the estimated values of T^s if the D values affected 
by the abundance He/H^ 0.2 are used in Tsff(D) relations cal- 
culated for He/H^ 0.1 at temperatures ranging form 23 000 K 
to 30 000 K. 

On average, due to the measurement errors of D, the ob- 
tained Tfiff(A[,D) parameters are obtained within uncertainties 
that for all studied early spectral types and luminosity classes 
are of the order of AT^slTes - 0.05, which is roughly the same 
as for the effective temperatures determined with the bolomet- 
ric flux method. Since the temperatures obtained with the bolo- 
metric flux method are not strongly sensitive to the character- 



istics of model atmospheres, they may be taken as a reference 
to study the properties of stellar atmospheres using theoretical 
spectroscopy. 

The {A\,D) parameters are less easily obtained for very hot 
stars. However, for active stars, like Be, B[e] and He-peculiar 
OB-types, the determination of BCD parameters might be as 
straightforward as it is for OB-type stars without emission. This 
is due to the fact that the photospheric component of BD is sep- 
arated from that originating in the circumstellar environment. 
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Appendix A: About the BCD system 



0.4 



The BCD parameters 
(Ai,D,$,b,^uv) 
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Fig.A.l. Graphical explanation of the BCD (Ai, D,<^rb,'^uv) 
pai^ameters. 

Initially, the BCD (Barbier-Chalonge-Divan) system has 
been thought as a two-parameter system of stellar classifica- 
tion, based on the spectrophotometric study of the continuum 
spectrum around the Balmer discontinuity. It is used for O, B, A 
and F-type stars. The two parameters are the size of the Balmer 
jump, D given in dex, and its mean spectral position, Ai given 
as the difference /li -3700 in A. 

The value of D is calculated at A = 3700 A, as D = 
logic ^3700+ /^3700- . where F3700+ is the Paschen side of the flux 



and F3700- is the flux in the Balmer continuum. The value of 
F3700+ is obtained by extrapolating the rectified Paschen con- 
tinuum to /t = 3700, for which a relation such as logFi/B^ = 
p X (1/A) + q is used. B,i can be the flux of a comparison star 
or simply the Planck function of a higher effective temperature 
than that expected for the studied star. In the present work we 
adopt the Planck function, so that the expression for D is; 



D = log 



37001 



If. 



3700 



/B 



dex 



(A.l) 



3700 



For the empirical determination of D, spectra of low reso- 
lution are used: AA - 8 A at A -3100. 

The average spectral position of the Balmer discontinuity is 
given by the point of intersection between the curve that passes 
over the maxima of fluxes in the limit of the Balmer line series, 
and the flux curve determined by the points log - D/2 on 
the Paschen side and log + D/2 on the Balmer side. In Fig. 
lA.ll we show the method of determining D and Ai. Since the 
wavelength scale to determine the intersection of flux curves 
is based on the intrinsic wavelengths of the identified spectral 
lines, the parameter /ii is not affected by any displacements due 
to the radial velocity of stars. 

When earring out the spectrophotometric study of the en- 
ergy distribution near the Balmer discontinuity, two other pa- 
rameters are obtained: the color gradient Ouv, given in fim and 
defined for the 3200-3700 A spectral region, and the Paschen 
gradient defined in two versions: <l>b or Oib for the spectral re- 
gions 4000-4800 A and 4000-6700 A, respective ly, both given 
in jum. A color gradient is defined as jAUe nl[l976l) : 



= 5^ 



d\nF] 



dd/A) 

which for a black body at temperature T, becomes : 
<l)(r) = (C2/r)(l-e-C2Mr^^ 



(A.2) 



(A.3) 



where C2 = hc/k - 1.4388 cm deg is the radiation constant. 
Assuming that for a given stellar energy distribution F^ it is 
O = const, between wavelengths Aa and Ah, the expression for 
O is: 



(A.4) 



The color gradient Ob was introduced in 1955 in the BCD 
system to distinguish F-type stars from B-type stars having the 
same {Ai,D) pairs. This parameter also has been used as a third 
BCD quantity related to the metal abundan ce of late type stars, 
in par ticular to the abundance ratio [Fe/H] ( Chalonge & DivarJ 
1977h . 

As the local temperature of the formation region of the 
Paschen continuum is close to the effective temperature, from 
( IA.4I ) and iA.3i we note that stars with the same effective tem- 
perature but different surface gravity defi ne a common region 
in the plane (/1 1 , D). This fact was used by lBarbier & Chalonge 
(119411) and Ichalonge & Divarj ( 1952 ) to determine the curvi- 
linear quadrilaterals that characterize the BCD classification 
system. To this end, the authors used the MK classification of 
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Fig.B.l. Error committed on the estimate of the observed frac- 
tion of the bolometric flux as a function of the uncertainty on 
the ISM color excess E{B - V). 



stars made by Morgan and Keenan d Keenan & MorganI 1951 ) 
themselves. They simply delimited the common region occu- 
pied by stars of the same MK spectral type with curves of in- 
trinsic constant Orb parameters. The same technique was used 
to draw the 'horizontal' lines that separate the MK luminosty 
classes. The BCD authors attempted to keep inside a common 
strip stars of all spectral types, but having the same MK lumi- 
nosity class label assigned by Morgan and Keenan. 

The color gradients Ob and Orb can be written as a func- 
tion of the {B — V) c olor of the UBV photometric system 
( lMouitahidetalJll998l) . The relation between the color ex- 
cesses i n the BCD and UBV syst ems due to the ISM reddening 
is then dChalonge & Divanlll973h : 



Av = 3.1£(B- V) = 1.7(0rb-0;?b) = 1.9(0b-0;;) mag , (A.5) 

where 0°j^ is the stellar intrinsic color gradient. 

One of the greatest advantages of the BCD method is that 
low resolution spectra are used. Since the BCD parameters can 
be easily obtained through automatic treatment of data, they 
are convenient for studies of stars in clusters or other stel- 
lar systems observed with multi-object spectrographs and/or 
spectro-imaging devices. The BCD system can also be of inter- 
est to analyze stars in far away galaxies that will be observed 
with the Extremely Large Telescopes. 

Appendix B: Calculation of the random errors 
affecting T^^ and 9^ 

As quoted in ^5.21 the T^^^ and 6^ parameters derived in the 
present work have the following sources of error: a) the ISM 
color excess E{B-V)\ b) the line blocking enhancement param- 
eter y; c) the log g parameter on which the model fluxes used to 
estimate the angular diameter depend and the filling factor due 
to the unobserved spectral region; d) the filling factor 5 used in 
relation (O to calculate the bolometric flux. 

a) The ISM color excess E{B - V). The correction of the ob- 
served energy distributions for the ISM extinction is probably 



the source of eiTor that has the heaviest consequences on the 
determination of T^^.. Both the value of the bolometric flux / 
in Eq. (O and the flux ratio f^ / entering Eq. depend on 
this correction. 

We write the error on the estimate of the color excess 
E(B - V) as AE = E(B - V) - E„{B - V), where E^iB - V) 
is the 'unknown' correct value of this excess. The effect of AE 
on the integrated fluxes coming from the observed spectral re- 
gion, /obs, can be estimated with model atmospheres, by calcu- 
lating the ratio A /obs//obs ^L/obsCAfi) -/obs]//obs- The behavior 
of A /obs //obs against AE for several eff'ective temperatures is 
shown in Fig. IB. II This calculation is made only for log g = 3.0 
since the effect of log g on A /obs /./obs is negligible. 

The uncertainty AE also affects the monochromatic fluxes 
f° used to derive the angular diameter. The angular diameter 
affects r^jj, which, in turn, determines T^siy) and the model 
flux entering relation For an estimate of the ratio AO^ jO^ 
we can simply use A - 0.1 pm, the middle wavelength of the 
interval over which 6^ is calculated. To make these calculations 
easier, the flux is represented as a function of T^siy) and 
log g, using the following interpolation expression: 



FMM[T,«(y),g] = A(logg)reff(r)*('°^«' 

A{\ogg) = 5.383 - 2.351 logg-H 0.32073 log^g (B.l) 
B(log^) = 1. 405 -i- 0.1 9 Hogg- 0.027 log^g. 

For the observed stellar flux at /I = 0.1 corrected for 
the ISM extinction is /°^ =/°^nO° ''*^»'^<*-^>], with A;o.7 =2.39, 
using ( IB. Il l the total uncertainty affecting the angular diameter 
can be written to the first order as: 



f = l.lAE{B-V)-\[B(\ogg)^ + [ 



lnreff(r)1gf^]Alogg; 



a In 4 (log g) 
fllogg 



(B.2) 



where it also includes the dependence on A log g and the uncer- 
tainty Sy in the factor A Tsff{y)/Tgff(y), whose dependence on 
6y we shall discuss below. 

In this work, the uncertainty inherent to E{B - V) was char- 
acterized with the Icr dispersion of the independent estimates 
of E(B - V) for a given star, i.e. AE(B -V)- ±o-e(b-v)- For the 
studied stars, it is on average crg(g_v) - 0.03 mag. In ( IB. 21 ) the 
direct dependence of A0^ on AE{B - V) is noted. However, 
it can be shown that the change of f^j with AE(B - V) is 
larger than that of Fqj induced by the coiTesponding AT^^^. 
This means that if AE{B - V) caiTies an increase of and 
an overestimation of Tj^^ through a greater bolometric flux, the 
corresponding increase of Fq j will compensate for the increase 
of /q 7, which maintains 9^ almost unchanged. 

b) The line blocking enhancement parameter y. To some 
degree, the estimate of y is also dependent on the ISM extinc- 
tion, whose effects in the far-UV can be confused with those 
produced by the line blanketing. The effect of the error on the 
enhancing factor y can be important as it enters the determina- 
tion of the temperature 7"eff(y) used to calculate both 9 and S. 
Nevertheless, the total effect of y on the value of T^^^ does not 
seem to carry larger deviations than AFetf ~ 300 K, as shown 
in Fig.|2] 
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The relative error Areff(7')/reff(y) as a function of the un- 
certainty on y can be estimated using the relation ( fT2] ). Since 
an interpolation expression of the function A as a function of 
T^eff (t) and log g is not easy to derive, to a good first order we 
can attempt an estimate of the uncertainty on T^sij) produced 
by Ay, with: 

- TA[reff(r), iogg]Ar . (B.3) 

The model fitting to the observed energy distributions car- 
ries uncertainties on the y of the order of Ay ^ 0.25. The con- 
sequences of this uncertainty on the value of Tttfiij) for a given 
r^jj. can be derived using ( IB. 31 ) or by interpolation in Table [T] 

c) The log g parameter. We assumed that the adopted model 
energy distributions correspond to gravities determined within 
an error A log g = 0.5 dex, which is twice the typical dispersion 
seen in Fig.|4] 

d) The filling factor 6. The bolometric flux fraction 6 would 
be realistic, if the actual stellar atmosphere of the observed 
star obeyed the basic assumptions that rule the models used. 
Although the non-LTE BW models produce somewhat differ- 
ent energy distributions in the far- and extreme-UV spectral 
regions, they cannot carry any significant change to the esti- 
mate of 6, as compared to the values derived with wind-free 
models. Having no definite knowledge of the phenomena that 
models might still be missing, we assumed that we could have 
uncertainties affecting 6 of the order of AS /S = 0.10. 

Since the uncertainties mentioned in this section do not 
propagate symmetrically (i.e. see Table lU, their combined ef- 
fect was simulated as follows. Let us call |Ap| the absolute error 
assigned to a given P quantity, i.e. P = E{B - V), y, logg and 
6. We divide each |Ap| into 4 parts: ep = |Ap|/4, so that we 
can assign to each quantity P 9 independent estimates, i.e. P = 
P - 4ep, P - 3ep,..., P,..., P + 3e/>, P + 4ep. Then, as we have 4 
different parameters to which we assign an error, each of them 
is taken with one possible simulated value P + nep (n from 1 
to 4); we have 9"^ = 6561 possible combinations that produce as 
many estimates of T^^ and 0-^ around the 'central' one where 
all Ap = 0.0. The frequency of T^^ and 8-^ values thus obtained 
have a bell-shaped distribution, not always symmetric. Let us 
call 2(fr,.ir) and Qiee) these distributions, where ej-.i^ = T^ff - T°^ 
and €0 = 9 - SP are the displacements from the respective T°^ 
and OP values for which all Ap = 0.0. The errors indicated in 
Table|2]are then the mean errors defined as: 

X+oo p+oo 
kie(f)de/ J e(e)de (B.4) 

that would correspond to e = yITJncr if the distributions Q{e) 
are 'normal', Gaussian, with dispersion cr. On average it is 
(?/rf^,) ^ 0.05, which is of the same order as the error expected 
for the Tfis(A\,D) values estimated in ^7.21 c.f. relation (l24l l. 
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Table 2: Program stars, observed and derived parameters. 



HD 


SpT/LC 




D 


E{B-V) 






r 










dex 


mag 


K 


mas 




K 


358 


B5.5V 


58.0 


0.288 


0.023 


12750±590 


1.016+0.010 


1.00 


14580 


886 


B2IV 


54.0 


0.155 


0.007 


21870±1390 


0.444+0.021 


1.00 


22740 


2905 


Blla 


31.0 


0.070 


0.339 


22160±1420 


0.348+0.021 


1.13 


21330 


3360 


B2IV 


54.0 


0.161 


0.031 


21850±1390 


0.307+0.021 


1.00 


22220 


3369 


B5V 


54.0 


0.259 


0.046 


15980+910 


0.305+0.015 


1.00 


16030 


4142 


B5.5V 


59.0 


0.283 


0.044 


16430±940 


0.162+0.015 


1.00 


14720 


4727 


B4V 


58.0 


0.249 


0.029 


16290±930 


0.266+0.015 


1.00 


16460 


11241 


B2V 


65.0 


0.162 


0.066 


22080±1410 


0.141+0.021 


1.00 


22420 


12767 


B5V 


55.1 


0.279 


0.000 


13250±640 


0.281+0.011 


1.00 


15000 


12953 


A2Ia 


-1.0 


0.231 


0.599 


9750±330 


0.591+0.007 


1.58 


9940 


13267 


B6Ia 


16.0 


0.155 


0.490 


14380±760 


0.259+0.012 


0.76 


13460 


14055 


AlV 


71.0 


0.471 


0.019 


9440±310 


0.540+0.007 


1.00 


9840 


14228 


B6V 


61.7 


0.296 


0.015 


12470±560 


0.507+0.010 


1.00 


14150 


14489 


AOIb 


11.0 


0.310 


0.420 


9840±340 


0.549+0.007 


2.04 


10650 


14818 


B2Ib-Ia 


31.0 


0.086 


0.496 


18300±1100 


0.202+0.018 


1.32 


20260 


15130 


A0.5IV 


52.8 


0.490 


0.019 


9920±350 


0.346+0.008 


1.00 


10040 


15318 


B9.5V 


61.9 


0.440 


0.020 


10630±400 


0.421+0.008 


1.00 


10700 


16046 


AOIV 


56.4 


0.448 


0.021 


10910±420 


0.314+0.008 


1.00 


10730 


16582 


B2IV 


52.6 


0.154 


0.020 


23160±1500 


0.239+0.022 


1.00 


22610 


16908 


B4V 


54.0 


0.231 


0.046 


17520±1030 


0.243+0.017 


1.00 


17430 


17081 


B6IV 


46.0 


0.311 


0.016 


13680±680 


0.344+0.011 


1.00 


14150 


17573 


B7V 


54.0 


0.339 


0.014 


12980±610 


0.476+0.010 


1.00 


12950 


18604 


B6V-IV 


50.4 


0.301 


0.033 


13940±710 


0.286+0.011 


1.00 


14380 


19356 


B7V 


60.0 


0.324 


0.047 


12800±590 


1.033+0.010 


1.00 


13240 


20041 


AOIa 


10.0 


0.276 


0.816 


10800±420 


0.643+0.008 


2.02 


10980 


21291 


B9Ia 


9.0 


0.213 


0.495 


11420±470 


0.794+0.009 


1.28 


11640 


21364 


B7.5V 


64.8 


0.337 


0.044 


13070±620 


0.473+0.010 


1.00 


12510 


21389 


AOIa 


4.0 


0.233 


0.633 


11040±440 


0.886+0.009 


1.22 


10860 


21447 


AlV 


72.0 


0.466 


0.026 


9240±300 


0.337+0.007 


1.00 


9890 


21790 


B9IV 


45.8 


0.408 


0.017 


11760±500 


0.317+0.009 


1.00 


11610 


21856 


BlIV 


57.0 


0.116 


0.183 


25370±1670 


0.124+0.023 


1.00 


27200 


22928 


B5III 


42.0 


0.281 


0.044 


14890±820 


0.583+0.013 


1.00 


15100 


22951 


BIV 


62.0 


0.112 


0.253 


29330±1980 


0.183+0.024 


1.00 


27990 


23180 


BlIII 


49.0 


0.121 


0.267 


22840±1470 


0.384+0.021 


0.97 


24190 


23227 


B4.5III 


43.4 


0.247 


0.019 


16230±930 


0.215+0.015 


1.00 


16470 


23288 


B7V 


56.0 


0.329 


0.117 


14020+720 


0.224+0.011 


1.00 


13220 


23324 


B7V 


56.0 


0.334 


0.061 


13210±640 


0.201+0.010 


1.00 


13050 


23338 


B61V 


46.0 


0.296 


0.059 


14180+740 


0.351+0.011 


1.00 


14650 


23408 


B6.5III 


38.5 


0.313 


0.129 


14310+750 


0.464+0.011 


1.00 


13960 


23625 


B2V 


59.0 


0.175 


0.294 


23980+1560 


0.114+0.022 


1.00 


21350 


23753 


B7IV 


50.0 


0.350 


0.065 


12680+580 


0.228+0.010 


1.00 


12740 


23850 


B7III 


39.0 


0.339 


0.086 


13020+620 


0.526+0.010 


1.00 


13240 


23923 


B9V 


58.0 


0.411 


0.054 


11640+490 


0.170+0.009 


1.00 


11350 


24131 


BIV 


61.0 


0.124 


0.265 


29130+1960 


0.132+0.024 


1.00 


26670 


24398 


Bllb 


37.0 


0.085 


0.328 


22040+1410 


0.659+0.021 


0.98 


22580 


24431 


09.5III 


54.0 


0.077 


0.653 


27370+1830 


0.155+0.023 


1.41 


31380 


24760 


BOIV-III 


54.5 


0.089 


0.086 


27160+1810 


0.394+0.023 


1.00 


29810 


25204 


B4IV 


50.1 


0.243 


0.050 


16970+990 


0.447+0.016 


1.00 


16829 


25490 


>A1V 


74.0 


0.496 


0.007 


8990+280 


0.605+0.007 


1.00 


9300 


27376 


B6V 


61.5 


0.303 


0.012 


12460+560 


0.522+0.010 


1.00 


13910 


27396 


B4IV 


51.0 


0.240 


0.146 


16720+970 


0.267+0.016 


1.00 


16960 


27962 


AlV 


71.0 


0.463 


0.000 


8680+260 


0.515+0.008 


1.00 


9970 
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et al.: Fundamental parameters of B Supergiants from the BCD System 
Table 2: continued. 



HD 

11 LJ 


op 1 / 


/I, 


LJ 


11 { I J V ) 


T a + AT a 

i eff ± eff 


u 31 ZAC/ 


y 




29248 


B2IV 


51.2 


0.144 


0.048 


23610±1530 


0.262+0.022 


1.00 


23120 


29305 


B7V 


61.3 


0.328 


0.021 


12120±530 


0.603+0.010 


1.00 


13030 


30211 


B4IV 


45.0 


0.248 


0.024 


15770±890 


0.346+0.015 


1.00 


16510 


30614 


B0.5Ia 


32.0 


0.057 


0.277 


25340±1670 


0.274+0.023 


0.95 


22620 


31647 


AlV 


76.0 


0.455 


0.032 


9750±330 


0.346+0.007 


1.00 


9950 


32309 


B9V 


63.0 


0.425 


0.018 


10800±420 


0.315+0.008 


1.00 


10960 


32549 


B9IV 


45.1 


0.438 


0.000 


9520±320 


0.392+0.007 


1.00 


11040 


32630 


B4V 


59.0 


0.232 


0.020 


17940±1070 


0.460+0.017 


1.00 


17390 


33802 


B6V 


59.5 


0.299 


0.028 


13280+640 


0.336+0.011 


1.00 


14140 


33904 


B7IV 


49.6 


0.331 


0.015 


12390±550 


0.590+0.010 


1.00 


13390 


33949 


B7III 


41.3 


0.358 


0.033 


12750±590 


0.358+0.010 


1.00 


12710 


34085 


B8Ia 


12.0 


0.177 


0.058 


12130±530 


2.713+0.010 


1.18 


12330 


34503 


B6IV 


44.8 


0.303 


0.045 


14450±770 


0.470+0.012 


1.00 


14440 


35468 


B2III 


49.0 


0.146 


0.023 


21840±1390 


0.775+0.021 


1.00 


22530 


35497 


B5.5IV 


44.0 


0.288 


0.022 


13960+710 


1.116+0.011 


1.00 


14910 


35600 


AOIb 


18.0 


0.375 


0.288 


10870±420 


0.324+0.008 


0.69 


10260 


36267 


B5V 


60.9 


0.257 


0.020 


16020+910 


0.314+0.015 


1.00 


15930 


36371 


B5Ia 


15.0 


0.129 


0.471 


15370+860 


0.476+0.014 


1.17 


14060 


36512 


08-9V 


68.1 


0.079 


0.029 


32340+2220 


0.142+0.025 


1.00 


32060 


36822 


BOIV 


55.0 


0.095 


0.111 


28340+1900 


0.198+0.024 


1.00 


29330 


37128 


BOIb 


39.1 


0.060 


0.059 


24670+1620 


0.681+0.022 


1.16 


24440 


37468 


09.5V 


71.0 


0.082 


0.061 


31270+2130 


0.230+0.025 


1.00 


31420 


37481 


B2V 


61.0 


0.142 


0.035 


24350+1590 


0.101+0.022 


1.00 


24700 


37744 


B2V 


65.0 


0.141 


0.057 


25690+1700 


0.088+0.023 


1.00 


25000 


38666 


08-9V 


66.5 


0.072 


0.017 


31510+2150 


0.111+0.025 


1.00 


33170 


38771 


BOII 


46.0 


0.081 


0.057 


23170+1500 


0.620+0.022 


1.00 


26730 


39970 


B9Ib-Ia 


12.0 


0.261 


0.487 


11530+480 


0.349+0.009 


1.51 


11250 


40111 


Bllb 


38.0 


0.084 


0.183 


24660+1620 


0.197+0.022 


1.22 


22950 


40183 


AlV 


69.0 


0.485 


0.012 


8910+280 


1.506+0.007 


1.00 


9680 


40312 


B9IV 


47.0 


0.448 


0.000 


9890+340 


0.959+0.008 


1.00 


10850 


40589 


B9Ib 


12.0 


0.284 


0.378 


11660+490 


0.295+0.009 


1.33 


10970 


41117 


B2Ia 


24.0 


0.050 


0.494 


20740+1310 


0.371+0.020 


1.53 


20640 


41753 


B4V 


55.0 


0.229 


0.032 


18800+1140 


0.251+0.018 


1.00 


17560 


42087 


B3Ia 


25.0 


0.102 


0.384 


16460+940 


0.258+0.015 


1.18 


17440 


43112 


BIV 


58.0 


0.117 


0.032 


27900+1870 


0.091+0.024 


1.00 


27190 


43384 


B4Ia 


20.0 


0.125 


0.621 


14780+800 


0.304+0.013 


1.38 


15140 


44743 


BlIII 


45.0 


0.112 


0.036 


25320+1670 


0.586+0.023 


1.00 


23690 


46300 


Allb 


17.0 


0.421 


0.083 


9800+340 


0.455+0.007 


1.00 


9740 


46769 


B7III-II 


30.0 


0.288 


0.151 


13930+710 


0.210+0.011 


0.45 


12990 


47105 


A1.5IV 


61.0 


0.518 


0.001 


9040+280 


1.435+0.007 


1.00 


9280 


47240 


Bllb 


35.0 


0.090 


0.372 


21540+1370 


0.157+0.021 


1.33 


21670 


47432 


08-9Ib 


43.0 


0.050 


0.389 


25620+1690 


0.133+0.023 


1.77 


28550 


47670 


B7III 


33.0 


0.331 


0.014 


12120+530 


0.625+0.010 


1.09 


12760 


47964 


B7IV-III 


42.0 


0.347 


0.021 


12100+530 


0.193+0.009 


1.00 


13040 


48434 


BOII 


44.0 


0.072 


0.241 


25290+1670 


0.128+0.023 


1.87 


26770 


48977 


B3V 


60.0 


0.210 


0.039 


19590+1210 


0.125+0.019 


1.00 


18820 


49567 


B4III 


37.0 


0.209 


0.051 


17270+1010 


0.125+0.016 


1.11 


17110 


52089 


BlIII 


45.0 


0.120 


0.034 


22010+1400 


0.801+0.021 


1.30 


23240 


53138 


B4Ia 


22.0 


0.130 


0.038 


14920+820 


0.580+0.013 


1.00 


15620 


53244 


B6III 


41.9 


0.289 


0.033 


13690+690 


0.378+0.011 


1.00 


14830 


58350 


B5Ia 


19.0 


0.152 


0.012 


12670+580 


0.882+0.010 


1.00 


14180 


66811 


06-7III 


55.3 


0.051 


0.043 


37250+2630 


0.371+0.027 


1.00 


36630 


67797 


B4IV 


50.9 


0.251 


0.028 


16680+960 


0.280+0.016 


1.00 


16450 


68520 


B5IV 


46.0 


0.287 


0.022 


14090+730 


0.326+0.011 


1.00 


14950 
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et al.: Fundamental parameters of B Supergiants from the BCD System 
Table 2: continued. 



OJ-/ 


op 1 / 


/I, 


LJ 


CyD V ) 


T a + AT a 

i eff ± ^-f eff 




y 




71155 


AlV 


74.0 


0.467 


0.025 


10060±360 


0.538+0.008 


1.00 


9800 


74280 


B2IV-V 


53.3 


0.189 


0.019 


19410±1200 


0.254+0.019 


1.00 


20080 


77327 


Aim 


47.0 


0.509 


0.008 


9080±290 


0.691+0.007 


1.00 


9790 


79447 


B3IV 


50.5 


0.220 


0.015 


17250±1010 


0.333+0.016 


1.00 


18010 


79469 


AOV 


63.7 


0.441 


0.023 


10980±430 


0.490+0.009 


1.00 


10630 


83754 


B5V 


59.1 


0.253 


0.024 


16150±920 


0.211+0.015 


1.00 


16220 


83944 


B9V 


66.2 


0.410 


0.006 


11550±480 


0.335+0.009 


1.00 


11150 


86440 


B6II 


25.0 




0.049 


13980±720 


0.493+0.011 


1.18 


13650 


87737 


AOIb 


18.0 


0.411 


0.053 


9820±340 


0.685+0.007 


1.00 


9900 


89021 


AlV 


64.0 


0.488 


0.006 


8790±270 


0.754+0.007 


1.00 


9730 


91316 


B1.5Ib 


34.8 


0.094 


0.056 


19830±1230 


0.312+0.020 


1.21 


21370 


95418 


>A1V 


70.0 


0.496 


0.004 


9470+310 


1.112+0.007 


1.00 


9480 


97633 


>A1V 


67.0 


0.510 


0.022 


9180±290 


0.772+0.007 


1.00 


9280 


98664 


B9IV 


54.4 


0.438 


0.023 


10680±410 


0.466+0.008 


1.00 


10970 


98718 


B5V 


60.1 


0.261 


0.014 


16760+970 


0.346+0.016 


1.00 


15730 


100600 


B3V 


58.0 


0.218 


0.033 


18230±1090 


0.129+0.018 


1.00 


18280 


100841 


Aim 


34.6 


0.518 


0.020 


9880±340 


0.768+0.008 


1.00 


9720 


100889 


B9V 


57.2 


0.435 


0.021 


11280+460 


0.334+0.009 


1.00 


10940 


106625 


B7IV 


49.4 


0.348 


0.009 


12360±550 


0.799+0.010 


1.00 


12820 


106911 


B4V 


58.2 


0.238 


0.026 


15330±850 


0.321+0.014 


1.00 


17060 


108767 


A0.5V 


71.0 


0.448 


0.020 


10580+400 


0.777+0.008 


1.00 


10270 


109026 


B4V 


58.9 


0.227 


0.027 


16740±970 


0.357+0.016 


1.00 


17700 


111123 


BlIV 


54.3 


0.125 


0.041 


27030+1800 


0.780+0.023 


1.00 


25740 


112185 


AlIV 


58.0 


0.511 


0.000 


9240+300 


1.504+0.007 


1.00 


9530 


112413 


B7V 


59.0 


0.335 


0.021 


11630+490 


0.718+0.009 


1.00 


12890 


116656 


AlV 


67.0 


0.470 


0.033 


9340+300 


1.361+0.007 


1.00 


9980 


120315 


B4V 


61.0 


0.227 


0.015 


17870+1060 


0.832+0.017 


1.00 


17710 


123299 


AOIV 


53.0 


0.482 


0.025 


10430+390 


0.571+0.008 


1.00 


10190 


129056 


B2IV 


47.8 


0.145 


0.041 


23100+1490 


0.537+0.022 


1.00 


22380 


129246 


>A1V 


79.0 


0.488 


0.000 


8990+280 


0.611+0.007 


1.00 


9220 


132058 


B2IV 


52.1 


0.148 


0.024 


24090+1570 


0.436+0.022 


1.00 


23000 


135742 


B7.5IV 


44.0 


0.363 


0.022 


12300+550 


0.801+0.010 


1.00 


12500 


137422 


>A2III 


35.0 


0.562 


0.022 


8280+240 


1.022+0.009 


0.48 


9110 


139006 


AlV 


69.0 


0.476 


0.026 


9900+340 


1.175+0.008 


1.00 


9820 


141003 


AlV 


69.0 


0.479 


0.018 


8810+270 


0.682+0.007 


1.00 


9770 


144217 


BIV 


65.5 


0.115 


0.210 


30540+2080 


0.505+0.024 


1.00 


27750 


145389 


B8.5V 


65.0 


0.363 


0.040 


11700+490 


0.407+0.009 


1.00 


12030 


147394 


B5V 


53.0 


0.262 


0.034 


16350+940 


0.356+0.015 


1.00 


15900 


148112 


AOIV 


57.0 


0.458 


0.072 


9810+340 


0.417+0.007 


1.00 


10520 


149438 


BOIV 


58.4 


0.095 


0.040 


31440+2150 


0.333+0.025 


1.00 


29640 


149881 


B0.5III-II 


45.0 


0.094 


0.070 


23420+1520 


0.063+0.022 


1.35 


24720 


155125 


>A1V 


75.0 


0.478 


0.000 


8620+260 


1.205+0.008 


1.00 


9570 


155763 


B6IV 


45.0 


0.314 


0.026 


13420+660 


0.586+0.011 


1.00 


14080 


158094 


B7V 


53.4 


0.343 


0.020 


12360+550 


0.508+0.010 


1.00 


12840 


159975 


B8III 


35.5 


0.360 


0.243 


12790+590 


0.428+0.010 


0.94 


12470 


160578 


B1.5III 


50.1 


0.131 


0.033 


24720+1620 


0.492+0.022 


1.12 


23780 


160762 


B3IV 


51.0 


0.211 


0.036 


19100+1170 


0.330+0.019 


1.00 


18590 


164353 


B5Ib 


23.8 


0.181 


0.185 


15420+860 


0.045+0.014 


1.11 


14790 


166182 


B2.5III 


43.0 


0.174 


0.074 


22420+1440 


0.232+0.021 


1.00 


19870 


169022 


AOIII 


45.5 


0.498 


0.024 


9520+320 


1.468+0.007 


1.00 


9960 


172167 


AlV 


66.0 


0.489 


0.000 


9470+310 


3.292+0.007 


1.00 


9670 


173300 


B7IV 


44.8 


0.332 


0.118 


14990+830 


0.597+0.014 


1.00 


13490 


175191 


B2.5V 


62.1 


0.187 


0.000 


18890+1150 


0.711+0.018 


1.00 


20370 


176437 


Aim 


32.0 


0.493 


0.040 


10000+350 


0.740+0.008 


1.00 


10000 
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et al.: Fundamental parameters of B Supergiants from the BCD System 
Table 2: continued. 



n LJ 


op 1 / 






£L\D V ) 


T a + AT a 

i eff ± eff 


ft + Afl 

U HI LxU 


y 




177724 


AlV 


66.0 


0.493 


0.052 


9830±340 


0.865+0.007 


1.00 


9600 


177756 


B9V 


58.9 


0.396 


0.004 


11780±500 


0.558+0.009 


1.00 


11600 


179761 


B8V-IV 


52.4 


0.364 


0.073 


13060±620 


0.268+0.010 


1.00 


12320 


182255 


B5V 


56.0 


0.268 


0.030 


15300±850 


0.208+0.014 


1.00 


15520 


186882 


AOIV 


50.0 


0.479 


0.031 


10150±360 


0.858+0.008 


1.00 


10280 


188209 


BOIb 


42.0 


0.065 


0.159 


25260±1670 


0.124+0.023 


1.33 


26040 


191692 


B9IV 


52.0 


0.462 


0.011 


10340±380 


0.712+0.008 


1.00 


10570 


192425 


AlV 


77.0 


0.461 


0.033 


9120±290 


0.370+0.007 


1.00 


9810 


192907 


AOV 


59.0 


0.451 


0.023 


10500±390 


0.408+0.008 


1.00 


10570 


195556 


B4II 


33.0 


0.204 


0.128 


17680±1040 


0.236+0.017 


1.07 


16320 


195810 


B6IV 


48.0 


0.303 


0.032 


14540±780 


0.370+0.012 


1.00 


14370 


196867 


B9IV 


53.0 


0.417 


0.018 


11220+450 


0.512+0.009 


1.00 


11370 


197345 


A2Ia 


4.0 


0.366 


0.056 


8720±260 


2.255+0.008 


1.70 


8630 


198001 


>A1V-IV 


64.0 


0.529 


0.021 


9370±310 


0.612+0.007 


1.00 


8980 


198478 


B4Ia 


20.0 


0.109 


0.562 


15390±860 


0.517+0.014 


1.21 


15880 


199081 


B4V-IV 


51.0 


0.248 


0.028 


16180±920 


0.240+0.015 


1.00 


16590 


202850 


B9.5Ia 


10.0 


0.265 


0.200 


11170±450 


0.536+0.009 


1.78 


11110 


204172 


BOII-Ib 


42.0 


0.073 


0.142 


24110±1570 


0.110+0.022 


1.36 


25370 


205021 


BIV 


59.0 


0.120 


0.037 


26920±1790 


0.312+0.023 


1.00 


26960 


206672 


B3ni 


43.0 


0.208 


0.089 


18360±1100 


0.241+0.018 


1.27 


18150 


207260 


A2Ia 


6.0 


0.378 


0.506 


8980±280 


1.010+0.007 


1.00 


8790 


207330 


B3III 


44.0 


0.190 


0.078 


17890±1060 


0.301+0.017 


1.00 


19270 


207971 


B7IV 


47.1 


0.338 


0.011 


12520±570 


0.649+0.010 


1.00 


13230 


209481 


06-7IV 


57.0 


0.062 


0.348 


30570±2080 


0.149+0.024 


1.00 


34280 


209744 


BlIII 


50.0 


0.120 


0.363 


25440±1680 


0.112+0.023 


1.00 


24510 


209819 


B8V 


64.3 


0.348 


0.020 


12310±550 


0.376+0.010 


1.00 


12320 


209952 


B5V 


51.7 


0.282 


0.013 


13920±710 


1.084+0.011 


1.00 


14990 


209961 


B2IV 


49.0 


0.180 


0.159 


21100±1330 


0.119+0.020 


1.00 


20420 


209975 


Blla 


32.0 


0.063 


0.331 


24720±1620 


0.222+0.022 


0.45 


22190 


210191 


B4III 


38.8 


0.202 


0.040 


17890±1060 


0.140+0.017 


1.08 


17750 


210418 


>A1V 


76.0 


0.480 


0.033 


8840±270 


0.747+0.007 


1.00 


9500 


212061 


AOV 


64.5 


0.458 


0.022 


10490±390 


0.527+0.008 


1.00 


10280 


212120 


B5V 


55.0 


0.282 


0.050 


15250±850 


0.288+0.014 


1.00 


14900 


212593 


B8.5II 


21.0 


0.299 


0.179 


11150±440 


0.441+0.009 


1.53 


11370 


212883 


B2.5IV-III 


48.0 


0.183 


0.085 


20400±1280 


0.101+0.020 


1.00 


20140 


212978 


B2III 


46.0 


0.176 


0.081 


20990±1320 


0.113+0.020 


1.00 


20360 


213420 


B2III 


44.0 


0.167 


0.131 


19750±1230 


0.267+0.019 


1.00 


20460 


213558 


A1.5V 


75.0 


0.465 


0.038 


9840±340 


0.590+0.007 


1.00 


9800 


213976 


B2V 


64.0 


0.165 


0.105 


22800±1470 


0.073+0.021 


1.00 


22170 


213998 


B9V 


58.5 


0.397 


0.017 


11740±500 


0.434+0.009 


1.00 


11590 


214240 


B5III 


43.0 


0.265 


0.126 


16330±930 


0.139+0.015 


1.00 


15750 


214652 


B2V 


60.0 


0.183 


0.105 


22350±1430 


0.081+0.021 


1.00 


20730 


214680 


06-7V 


65.0 


0.068 


0.089 


32380+2220 


0.136+0.025 


1.00 


33820 


214923 


B9IV 


50.3 


0.420 


0.005 


11430±470 


0.581+0.009 


1.00 


11360 


214993 


BlIII 


49.0 


0.120 


0.111 


24130+1570 


0.155+0.022 


1.00 


24260 


214994 


A0.5V 


62.0 


0.477 


0.050 


9930+350 


0.364+0.008 


1.00 


9970 


215191 


B2IV 


52.0 


0.145 


0.131 


22800+1470 


0.099+0.021 


1.00 


23230 


217101 


B2V 


63.0 


0.151 


0.093 


22760+1460 


0.107+0.021 


1.00 


23700 


217811 


B4m 


45.0 


0.225 


0.217 


19330+1190 


0.132+0.019 


1.00 


17420 


218045 


AlIV 


55.0 


0.500 


0.021 


9850+340 


1.036+0.008 


1.00 


9850 


218376 


B0.5III 


46.0 


0.095 


0.244 


27100+1810 


0.204+0.023 


1.28 


25190 


218407 


B3V 


60.0 


0.184 


0.174 


21200+1340 


0.101+0.020 


1.00 


20650 


219688 


B5V 


51.8 


0.266 


0.021 


15260+850 


0.297+0.014 


1.00 


15730 


222173 


B8IV 


44.0 


0.371 


0.038 


12620+580 


0.369+0.010 


1.00 


12340 



J. Zorec et al.: Fundamental parameters of B Supergiants from the BCD System 
Table 2: continued. 



HD SpT/LC 



D 



E(B - V) Teff + A^eff 



e + AG 



7 



222661 
223640 

224990 



B9.5V 
B6.5V 

B5V 



66.0 0.428 
55.9 0.321 

54.2 0.259 



0.021 
0.000 

0.032 



10860±420 
11740±500 

16100+910 



0.381±0.008 
0.260+0.009 

0.215+0.015 



1.00 
1.00 

1 .00 



10820 
13510 

16030 



The parameter Ai is given as Ai - 3700 A. 

The notation >SpT/LC means that the spectral type is cooler than the indicated one. 



J. Zorec et al.: Fundamental parameters of B Supergiants from the BCD System 

Table 5: Compilation of effective temperatures determined for the pro- 
gram stars by other authors. 



HD 


^ eff 


''eff 


''eff 


''eff 


^eff 


''eff 


''eff 


358 


12750 














886 


21870 




21990 






21250 




2905 


22160 






21600 




24000 


22500 


3360 


21850 




22210 










3369 


15980 


15520 












4142 


16430 


16670 












4727 


16290 


15850 












11241 


22080 


25350 












12767 


13250 














12953 


9750 














13267 


14380 










15000 




14055 


9440 








9490 






14228 


12470 


12900 






13080 






14489 


9840 














14818 


18300 










20000 


18250 


15130 


9920 








9890 






15318 


10630 








10600 






16046 


10910 














16582 


23160 




23570 










16908 


17520 


16850 






16940 






17081 


13680 


13210 


13040 




13170 


13100 




17573 


12980 


12570 






12400 






18604 


13940 


13240 












19356 


12800 


14290 












20041 


10800 














21291 


11420 










11500 




21364 


13040 


12650 






12540 






21389 


11040 














21447 


9240 














21790 


11760 














21856 


25370 


25880 












22928 


14890 




13940 




14200 






22951 


29330 


28970 












23180 


22840 


22860 












23227 


16230 














23288 


14020 














23324 


13220 


12760 












23338 


14180 














23408 


14310 














23625 


23980 


25530 












23753 


12680 














23850 


13020 














23923 


11640 














24131 


29130 


28380 












24398 


22040 




19870 


19900 




23000 




24431 


27370 














24760 


27160 


25760 












25204 


16970 














25490 


8990 








9260 


9250 




27376 


12460 














27396 


16720 




16800 










27962 


8680 










9000 





J. Zorec 



et al.: Fundamental parameters of B Supergiants from the BCD System 
Table 5: continued. 



tiu 


^eff 




■r(2) 


■r(3) 


t(4) 


t{5) 


t(6) 
^eff 


29248 


23610 


— 


23500 


21600 


_ 


_ 


_ 


29305 


12120 


— 


_ 


_ 


_ 


_ 


_ 


30211 


15770 


15540 


14950 




15630 


_ 




30614 


25400 


_ 


25000 


26300 


_ 


32500 


29000 


31647 


9750 


_ 


_ 


_ 


_ 


_ 


_ 


32309 


10800 


_ 


10200 


_ 


_ 


_ 


_ 


32549 


9520 


_ 


10200 


_ 


_ 


_ 


_ 


32630 


17940 


17020 


17580 


_ 


17050 


16380 


_ 


33802 


13280 


_ 


11860 


_ 


_ 


_ 


_ 


33904 


12390 


_ 


_ 


_ 


_ 


_ 


_ 


33949 


12750 


11690 


_ 


_ 


_ 


_ 


_ 


34085 


12130 


_ 


11660 


_ 


_ 


13000 


_ 


34503 


14450 


13510 


14260 


_ 


14620 


_ 


_ 


35468 


21840 


21040 


21760 


20900 


21940 


_ 


_ 


35497 


13960 


13650 


13820 


_ 


13490 


_ 


_ 


35600 


10870 


_ 


_ 


_ 


_ 


11000 


_ 


36267 


16020 


17100 


_ 


_ 


_ 


_ 


_ 


36371 


15370 


_ 




_ 


_ 


16500 


_ 


36512 


32340 


_ 


34350 


_ 


_ 


_ 


_ 


36822 


28340 


_ 


_ 


_ 


31120 


_ 


_ 


37128 


24670 


_ 


24960 


26400 


_ 


28500 


27250 


37468 


31270 


_ 


31560 


_ 


_ 


_ 


_ 


37481 


24350 


24550 


_ 


_ 


_ 


_ 


_ 


37744 


25690 


24830 


_ 


_ 


_ 


_ 


_ 


38666 


31510 


_ 


_ 


_ 


_ 


_ 


_ 


38771 


23170 


_ 


26390 


23100 


_ 


27500 


26250 


39970 


11530 


_ 


_ 


_ 


_ 


_ 


_ 


40111 


24660 


_ 


_ 


21400 


_ 


_ 


_ 


40183 


8910 


_ 


_ 


_ 


_ 


_ 


_ 


40312 


9890 


_ 


_ 


_ 


_ 


_ 


_ 


40589 


11660 


_ 


_ 


_ 


_ 


_ 


_ 


41117 


20740 


_ 


17460 


_ 


_ 


19500 


19000 


41753 


18800 


16970 


17410 


_ 


17640 


_ 


_ 


42087 


16460 


_ 


_ 


_ 


_ 


20500 


18000 


43112 


27900 


28510 


_ 


_ 


_ 


_ 


_ 


43384 


14780 


_ 


_ 


_ 


_ 


_ 


_ 


44743 


25320 


_ 


25500 


24700 


_ 


_ 


_ 


46300 


9800 


_ 


_ 


_ 


_ 


_ 


_ 


46769 


13930 


14520 




_ 








_ 


47105 


9040 


_ 


9260 


_ 


9380 


9150 


_ 


47240 


21540 


20610 


_ 


_ 


_ 


_ 


_ 


47432 


25620 


_ 




_ 


_ 


_ 


_ 


47670 


12120 


11910 


11610 


_ 


11950 


_ 


_ 


47964 


12100 














48434 


25290 


24270 












48977 


19590 


19770 












49567 


17270 


17100 












52089 


22010 




21830 


21000 








53138 


14920 




14760 






18500 


16000 


53244 


13690 




13600 










58350 


12670 




13150 






16000 


15000 


66811 


37250 




33760 


35100 








67797 


16680 




15460 











J. Zorec 



et al.: Fundamental parameters of B Supergiants from the BCD System 
Table 5: continued. 



tiu 


^eff 


^efF 


t(2) 


t(3) 


t(4) 


t{5) 


t(6) 
^eff 


68520 


14090 


14420 


_ 


_ 


_ 


_ 


_ 


71155 


10060 


_ 


_ 


_ 


_ 


_ 


_ 


74280 


19410 


18680 


18790 


_ 


18650 


_ 


_ 


77327 


9080 


_ 


_ 


_ 


_ 


_ 


_ 


79447 


17250 


17220 


_ 


_ 


_ 


_ 


_ 


79469 


10980 


_ 


_ 


_ 


10640 


_ 


_ 


83754 


16140 




15400 


_ 


_ 


_ 


_ 


83944 


11550 


10740 


_ 


_ 


_ 


_ 


_ 


86440 


13980 


13300 


_ 


_ 


_ 


_ 


_ 


87737 


9820 


_ 


_ 


_ 


_ 


_ 


_ 


89021 


8790 


_ 


_ 


_ 


9160 


9000 


_ 


91316 


19830 


_ 




20300 


_ 


_ 


22000 


95418 


9470 


_ 


9170 


_ 


9650 


9600 


_ 


97633 


9180 


_ 


_ 


_ 


9130 


9250 


_ 


98664 


10680 


_ 


_ 


_ 


_ 


_ 


_ 


98718 


16760 


_ 


_ 


_ 


_ 


_ 


_ 


100600 


18230 


19320 


_ 


_ 


_ 


_ 


_ 


100841 


9880 


_ 


10470 


_ 


_ 


_ 


_ 


100889 


11280 


_ 


_ 


_ 


10810 


_ 


_ 


106625 


12360 


_ 


12450 


_ 


_ 


_ 


_ 


106911 


15330 


_ 


_ 


_ 


_ 


_ 


_ 


108767 


10580 


_ 


_ 


_ 


_ 


_ 


_ 


109026 


16740 


_ 


15510 


_ 


_ 


_ 


_ 


111123 


27030 


_ 


27600 


_ 


_ 


_ 


_ 


112185 


9240 


_ 


_ 


_ 


_ 


_ 


_ 


112413 


11630 


_ 


_ 


_ 


_ 


_ 


_ 


116656 


9340 


_ 




_ 


_ 


_ 


_ 


120315 


17870 


17360 


16770 


_ 


17100 


16900 


_ 


123299 


10430 


_ 


_ 


_ 


_ 


10000 


_ 


129056 


23100 


20320 


_ 


_ 


_ 


_ 


_ 


129246 


8990 


_ 


_ 


_ 


_ 


_ 


_ 


132058 


24090 


_ 


23780 


22800 


_ 


_ 


_ 


135742 


12300 


12060 


_ 


_ 


12310 


12130 


_ 


137422 


8280 


_ 


_ 


_ 


_ 


_ 


_ 


139006 


9900 


_ 


_ 


_ 


_ 


_ 


_ 


141003 


8810 


_ 


_ 


_ 


_ 


_ 


_ 


144217 


30540 


_ 


_ 


_ 


_ 


_ 


_ 


145389 


11700 


_ 


10860 


_ 


_ 


_ 


_ 


147394 


16350 


15030 


15210 


_ 


15050 


15000 


_ 


148112 


9810 


_ 


_ 


_ 


_ 


_ 


_ 


149438 


31440 


_ 


_ 


31800 


29920 


_ 


_ 


149881 


23420 


21130 


_ 


_ 


29280 


_ 


_ 


155125 


8620 


_ 


_ 


_ 


_ 


_ 


_ 


155763 


13420 




12960 




13090 






158094 


12360 


11480 


12010 










159975 


12790 




11650 










160578 


24720 




25780 


24400 








160762 


19100 


17640 


17810 




17460 






164353 


15420 


13490 








16500 


15500 


166182 


22420 


21380 


20320 










169022 


9520 




9460 




9540 






172167 


9470 




9660 




9620 






173300 


14990 


12390 






12110 







J. Zorec 



et al.: Fundamental parameters of B Supergiants from the BCD System 
Table 5: continued. 







■r(l) 


t(2) 


t(3) 


t(4) 
^eff 


t{5) 


t(6) 
^eff 


175191 


18890 


_ 


18990 


_ 


_ 


_ 


_ 


176437 


10000 


_ 


_ 


_ 


9830 


9550 


_ 


177724 


9830 


_ 


9480 


_ 


9550 


9500 


_ 


177756 


11780 


11560 


11410 


_ 


11430 


_ 


_ 


179761 


13060 


_ 


_ 


_ 


_ 


13000 


_ 


182255 


15300 


_ 




_ 


_ 


_ 


_ 


186882 


10150 


_ 


9880 


_ 


10180 


_ 


_ 


188209 


25260 




31080 


27700 


_ 


_ 


_ 


191692 


10340 


_ 


_ 


_ 


10170 


_ 


_ 


192425 


9120 


_ 


_ 


_ 


_ 


_ 


_ 


192907 


10500 


_ 


10650 


_ 


_ 


10250 


_ 


195556 


17680 


_ 


18140 


_ 


_ 


_ 


_ 


195810 


14540 


13550 


13610 


_ 




_ 


_ 


196867 


11220 


_ 


10820 


_ 


10910 


_ 


_ 


197345 


8720 


_ 


_ 


_ 


_ 


_ 


_ 


198001 


9370 


_ 


_ 


_ 


9180 


9200 




198478 


15390 


_ 


14270 


_ 


_ 


18000 


17170 


199081 


16180 


16560 


_ 


_ 


_ 


_ 


_ 


202850 


11170 


_ 


_ 


_ 


_ 


_ 


11000 


204172 


24110 


_ 


_ 


25200 


_ 


28500 


28500 


205021 


26920 


_ 


25780 


26200 


_ 


_ 


_ 


206672 


18360 


17390 


17890 


_ 


17110 


_ 


_ 


207260 


8980 


_ 


_ 


_ 


_ 


_ 


_ 


207330 


17890 


17100 




_ 


_ 


_ 


_ 


207971 


12520 


12170 


11800 


_ 


11950 


_ 


_ 


209481 


30570 


_ 


_ 


_ 


_ 


_ 


_ 


209744 


25440 


_ 




_ 


_ 


_ 


_ 


209819 


12310 


12470 


11280 


_ 


_ 


_ 


_ 


209952 


13920 


13450 


13660 


_ 


13930 


_ 


_ 


209961 


21100 


21430 


_ 


_ 


_ 


_ 


_ 


209975 


24720 


_ 


_ 


_ 


_ 


32500 


_ 


210191 


17890 


_ 


_ 


_ 


_ 


_ 


_ 


210418 


8840 


_ 


_ 


_ 


_ 


_ 


_ 


212061 


10490 


_ 


_ 


_ 


10290 


_ 


_ 


212120 


15250 


14580 


_ 


_ 


14690 


_ 


_ 


212593 


11150 




9930 


_ 


_ 


_ 


11800 


212883 


20400 


22700 


_ 


_ 


_ 


_ 


_ 


212978 


20990 


20370 


_ 


_ 


20790 


_ 


_ 


213420 


19750 


19850 


20840 


_ 


20490 


_ 


_ 


213558 


9840 


_ 


_ 


_ 


9600 


_ 


_ 


213976 


22800 


_ 


_ 


_ 


_ 


_ 


_ 


213998 


11740 


11490 


11220 


_ 


11530 


_ 


_ 


214240 


16330 


17380 


_ 


_ 


_ 


_ 


_ 


214652 


22350 


21930 












214680 


32380 




33560 




33110 


33000 




214923 


11430 


11320 






11500 






214993 


24130 


23600 












214994 


9930 








9560 


9530 




215191 


22800 


22930 






24050 






217101 


22760 


23120 












217811 


19330 


19100 












218045 


9850 








9850 






218376 


27100 




27470 






29000 





J. Zorec et al.: Fundamental parameters of B Supergiants from the BCD System 
Table 5: continued. 



HD 


eff 


7^(1) 
-"eff 


7^(2) 
eff 


7^(3) 
eff 


7^(4) 
eff 


7^(5) 
eff 


7^(6) 
eff 


218407 


21200 


23230 












219688 


15260 




15210 




14910 






222173 


12610 


11850 


11870 




11870 






222661 


10860 














223640 


11740 














224990 


16100 








15330 







ri/ : Gulati et al. (1989),Castelli (1991) 

Tf^ : Underhill et al. (1979), Code et al. (1976), Malagnini et al. (1986) 
rf^ : Remie & Lamers (1981) 
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